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Abstract: This paper is a continuation of |ArXiv:0707.13"2^ where improved holo- 
graphic theories for QCD were set up and explored. Here, the IR confining geometries 
are classified and analyzed. They all end in a "good" (repulsive) singularity in the 
IR. The glueball spectra are gapped and discrete, and they favorably compare to 
the lattice data. Quite generally, confinement and discrete spectra imply each other. 
Asymptotically hnear glueball masses can also be achieved. Asymptotic mass ratios 
of various glueballs with different spin also turn out to be universal. Meson dynamics 
is implemented via space filling D4 — D4 brane pairs. The associated tachyon dy- 
namics is analyzed and chiral symmetry breaking is shown. The dynamics of the RR 
axion is analyzed, and the non-perturbative running of the QCD ^-angle is obtained. 
It is shown to always vanish in the IR. 
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1. Introduction and summary 

This paper is a direct sequel of , and the reader is guided there for a comprehensive 
introduction and summary of results of both papers. In the first part of this work, 
two of the authors establish and motivate a general 5D holographic setup to describe 
4D gauge theories with a large number of colors (large Nc). The setup described 
there constitutes a bottom-up approach, motivated in part from known features of 
5D non-critical string theory and in part by what we expect from QCD. 

The pure gauge dynamics is encoded holographically in the solution of a two- 
derivative action for the expected 5D fields: the 5D metric (dual to the YM stress 
tensor), a scalar (the dilaton, dual to Tr[F^]), and a pseudoscalar (the axion, dual 
to Tr[F A F]). The dilaton potential is expected to be non-trivial and is expected to 
obtain non-trivial contributions from the non-propagating four-form. In practice, the 
potential (and the associated superpotential) is in one-to-one correspondence with 
the QCD /3-function and is chosen in such a way as to reproduce known features 
(e.g. UV asymptotic freedom and IR confinement) of the gauge theory. This is what 
makes our approach phenomenological. 

The resulting backgrounds present an improvement over pre-existing models of 
"phenomenological holography", e.g. p|, among other advantages, the back- 
grounds we present incorporate the running of the coupling and asymptotic freedom; 
establish a one-to-one correspondence between the 5D geometry and the gauge the- 
ory parameters, namely the /9-function /3(A), and the dynamically generated IR scale 
Aqcd and do not require specifying the boundary conditions in the IR. Moreover, 
they provide a natural environment to study non-perturbative dynamical phenomena 
such as confinement, generation of the mass gap and chiral- symmetry breaking. 

Part I is devoted to establish and motivate the setup, and to analyze the per- 
turbative UV regime of the correspondence. The present work, on the other hand, 
focuses on the analysis of the non-perturbative regime. As one of our main results, we 
establish a relation between color confinement (i.e. an area law for the Wilson loop) 
and the properties of the geometry in the IR, and we show that confining backgrounds 
always exhibit a mass gap and generically a discrete spectrum. This is a nontrivial 
statement, as in our models there is no IR boundary (which would automatically 
guarantee both confinement and a mass gap). In most of this work we focus on the 
pure Yang-Mills sector, which we describe holographically by a 5D Einstein-Dilaton 
system. We discuss the addition of Nf matter flavors in the quenched approxima- 
tion Nf -C Nc, so that we can neglect the backreaction of the 5D fields dual to the 
operators containing quarks. 

The structure of the present paper is as follows. In Section |^ we give an overview 
of the setup discussed in We recall how asymptotic freedom demands the UV ge- 
ometry to be asymptotically AdS^ with logarithmic corrections. We then review the 
holographic dictionary, mapping field theory quantities to their geometrical counter- 
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parts. In particular, there is a one-to-one correspondence between the superpotential 
associated to the geometry and the exact field theory /3-function. In the last sub- 
section we analyze the IR behavior of space-times that have AdSs UV asymptotics, 
showing that backgrounds that are not conformal in the IR necessarily exhibit an IR 
singularity. 

In Section ^ we provide a complete characterization of 5D asymptotically AdS^ 
backgrounds that exhibit confinement in the IR. Here, as a definition of confinement, 
we require that the Wilson loop exhibits an area law behavior. We compute the Wil- 
son loop holographically, using the prescription of @], as the action of a classical 
string world-sheet with fixed UV boundary conditions. We show that confinement 
requires the scale factor to vanish sufficiently fast in the IR. We formulate this re- 
quirement in terms of the superpotential and the 4D /3-function. We then discuss the 
holographic computation of the 't Hooft loop, relevant for the potential between two 
color-magnetic charges, in order to investigate whether screening behavior is present. 

In Section 4 we perform a holographic computation of the (regularized) QCD 
vacuum energy. 

In Section ^ we discuss the dynamics of the 5D axion. This field does not 
backreact on the geometry in the large limit, however its normal modes give the 
spectrum of 4D pseudoscalar glueballs. In addition, its UV source is associated to 
the t/ieta-angle of YM and its IR properties are relevant for the discussion of the 
effective QCD ^^-parameter. We argue that this field must vanish in the IR. In this 
way a pseudo-scalar glueball condensate is dynamically generated. This suggests 
that the effects of the 6'-parameter are screened in the IR. 

In Section ^ we discuss the general features of the low-energy particle spectrum in 
our model, obtained from the fluctuations of the bulk flelds around the background. 
For scalar and tensor glueballs, in all conflning potentials the spectrum is gapped 
and discrete. With the exception of a certain class of "pathological" geometries 
in which the singularity is not screened, the particle spectrum can be computed 
unambiguously imposing normalizability of the eigenfunctions. We flnd potentials 
where the glueballs have a linear asymptotic spectrum, i.e. ~ n. We discuss 
the addition of flavor branes along the lines of ||^, where it was proposed that chiral 
symmetry breaking is correctly described by open string tachyon condensation. We 
show that this idea can be naturally implemented in our setup. We discuss the 
asymptotics of the tachyon background and the qualitative features of the spectrum 
of mesons. 

In Section 7 we discuss the parameters of the gravity backgrounds and their re- 
lations with the gauge theory parameters. We show that, once conflning asymptotics 
in the IR are imposed, the 5D backgrounds are completely specifled by the dilaton 
potential plus the arbitrary choice of a single integration constant, that only affects 
the overall energy scale. It is in one-to-one correspondence to a choice of Aqqd in 
the dual fleld theory. 
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In Section ^ we give some concrete examples, in which we specify the exact 
/3-function, solve numerically for the corresponding geometry and compute, again 
numerically, the glueball spectrum. We compare our results with the available lattice 
data, and with similar computations in other models. In particular, we show that the 
"linear confinement" background fits particularly well the lattice data. On the other 
hand, in hard- wall models like 0, or generically in models with a "quadratic" 
mass spectrum, the agreement is not as good. 

Some of the technical details can be found in the appendices. In particular, 
the reader interested in the details of the characterization of confining backgrounds, 
including their various geometric properties, is referred to Appendix 0. 

2. Building blocks of holographic QCD 

In this section we review some properties of the 5d backgrounds and their holographic 
interpretation. Some of these points where extensively discussed in [^. 

2.1 The 5D backgrounds 

As detailed in we take a "minimal" non-critical approach to holographic large Nc 
QCD type theories, in which the 5D string theory dual has, as low energy excitations, 
the duals of the lowest- dimension gauge invariant operators. In the pure glue sector 
these are: the five-dimensional metric g^u (dual to the YM stress tensor); a scalar 
field which we call the dilaton (dual to the YM operator TrF^) and an axion, dual 
to TtFF. We may ignore the axion when searching for the QCD vacuum solution 
as its contribution is subleading in l/Nc ||2^. It can be included in the sequel (this 
is discussed in section ^ as it does affect some important physics, in particular that 
of flavor singlet mesons. The scalar encodes the running of the YM coupling, and 
it is naturally identified with the 5D string dilaton. 

We should emphasize that we think of the 5D bulk theory as a (non-critical) 
string theory, not just gravity. However, we restrict ourselves to the two-derivative 
effective action, including a general dilaton potential, that contains also a subclass 
of higher a' terms as argued in . 

Therefore, the string-frame action describing the low-lying excitations is: 

Ss = M^Nl j d'xV^^ 

where we have introduced the 't Hooft coupling 

A = N,e^ . (2.2) 

It is related to the 't Hooft coupling of the gauge theory, up to a multiplicative 
constant. Vs{X) is an "effective potential" originating from integrating-out the non- 
dynamical 4-form, and other higher curvature corrections. We do not attempt 



R 
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here to derive Vs{X) from first principles. We determine certain of its properties by 
requiring that the geometry that follows from Vs reproduces some known features of 
the Yang-Mills dynamics. In particular, requiring UV asymptotic freedom constraints 
the asymptotics of the potential for small values of A. Requiring an area law for the 
Wilson loop on the other hand constrains the asymptotics of the potential for large 
A. From now on we also define for convenience a renormalized dilaton $ as 

A = e* . (2.3) 

We mostly work with the Einstein frame metric, 

9,. = e-l''{gs)^u, (2.4) 

for which the action reads: 



R - ^g'^-'d^^d,^ + r($) 



V{^) = e^*/V5($) . (2.5) 



In the large limit we assume, g^u and $ are independent of Nc- 
We search backgrounds of the form: 

g^^ = du^ + e^^'^^^ri.jdx'dx^ = e^^^'^ {dr^ + rjijdx'dx^) , $ = $(u), (2.6) 

where are the 4D space-time coordinates, and rjij = diag{—,+,+,+). We write 
the metric in two different coordinate systems, related by: 

^ = e-^("). (2.7) 
du 

We name the first set of coordinate system the domain-wall coordinates. The second 
set involving r will be called the conformal coordinates as the metric is explicitly 
conformally flat in this coordinate system. Throughout the paper, we will use a 
prime for d/du and a dot for d/dr. 

The independent Einstein's equations take the following form (in domain-wall 
coordinates) : 

$'2(«) = ~A"{u), V($) = 3A"{u) + 12A'^{u). (2.8) 
These equations can be written in first-order form in terms of a superpotential ($) 



dW 4 
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Given any scale factor A{u) such that A"{u) < 0, one can invert the relation be- 
tween $ and u using the first equation (|2.8|) and find a superpotential = 
— 4/9 This determines a potential, such that the given A{u) is a solution. 

This useful property allows to investigate the backgrounds under consideration start- 
ing directly from a parametrization of the metric, rather than the dilaton potential. 
In conformal coordinates, Einstein's equations ( |2.8| ) read: 

^^(r) = -^(^A{r)-A\r)y V{^) = e-^"^^'^ (3A{r) + 9A\r)y (2.11) 
or, in terms of the superpotential: 

$ = ^e^, A = --We^. (2.12) 
(i$ 9 

As shown in [|1|, asymptotic freedom in the gauge theory translates into an 
asymptotic AdS^ region^, where the dilaton asymptotes to — oo, and W Wq > 0: 

r u 1 

A{u) u/i + 0{\ogu) , $(«) log ---log(M) +0(1), M^-oo, 

(2.13) 

or, in conformal coordinates: 

v4(r) ~ -log r/£ + (-^1 , $ ~ - log(- log(rA)) + 0(1), r ~ £e"/^ ^ 0. 

VogrJ 

(2.14) 

The AdS scale i is fixed by the value of W{X) for A = 0: i = {9/4)Wo^; A is an 
integration constant whose meaning will be clarified in section |^. The subleading 
terms are also fixed, order by order by matching the /3-function coefficients. This 
corresponds to a dilaton potential of the form V"($) ~ Vq + Vie* + V2e^* + 0(e^*). 
Since A ~ e* is small in this region, we expect this potential to be generated by 
the full resummation of the a' expansion, and is to be interpreted as an "effective 
potential." The information of its weak coupling expansion maps on the field theory 
side, to the perturbative /5-function. 

In this paper we are mostly concerned with the opposite regime, i.e. the large 
A IR region. One of our goals is to find which solutions to (|2.8|), satisfying the 
UV asymptotics (|2.13|), provide an area- law for the Wilson loop, and what kind of 



potentials are necessary to realize those solutions. Before addressing this problem, we 
give some preliminary discussion about the holographic dictionary and the infra-red 
properties of 5D geometries. 



^We should note that "asymptotically AdSs" here has a different meaning that the similar term 
in the mathematics literature (see for example 0| and the references therein) . Here the corrections 
to the AdSs metric are also logarithmic, while there they are powers of the conformal coordinate r. 
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2.2 Holographic dictionary 

In order to exploit the gauge- gravity duality, we must first establish a dictionary 
between the 5D and 4D quantities. In particular we must identify the Yang Mills 
coupling and energy scale on the gravity side. For this we use the dictionary estab- 
lished in 

At a given position in the fifth dimension, the four-dimensional energy scale is 
set by the scale factor e^{u), as seen from eq ( [^.6| ). This, we argue, leads to the 
identification: 

logE^A{u). (2.15) 

Close to the AdS^ boundary, this reduces to the familiar identification E = 1/r. The 
correspondence ( p.l5| ) does not fix the absolute units of the 4D energy scale with 
respect to the scale factor. This is consistent with the observation that a constant 
shift in A{u) leaves Einstein's equations ( [2.8| ) invariant, and can be used to change 
the unit energy in a given background. 

Notice that the overall scale factor in (|2.15|) is in the Einstein frame. In pure 
AdS^ with a constant dilaton this distinction does not matter, but in our backgrounds 
the dilaton does not asymptote to a constant in the UV, therefore this clarification 
is needed. In particular, the Einstein frame scale factor has the important property 
of being monotonically decreasing with u (see Section 2.3). This property is not 
shared by the string frame metric. Monotonicity is crucial if we want our geometry 
to be dual to a single RG flow from the UV to the IR (and not, e.g, to two different 
UV theories that flow to the same IR). For a discussion on the string corrections to 



see 



The /9-function 

We identify^ the 4D 't Hooft coupling A = Qym^c as^ 

A = e*. (2.16) 



With the identification (|2.15| ), it follows that the /3-function of the 't Hooft 



coupling is related to 5D fields as: 



or, in terms of the phase-space variable X, introduced in 



X = — , /3 = SAX (2.18) 



^There are some ambiguities in this identification that are discussed in [Q. 



^As stated in 2.2, the string dilaton, = log^g, is related to $ by ^ = $ — logA^c- ^ is the 



appropriate variable to use in the large N,, limit. 
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The above definition is independent of reparametrizations of tlie radial coordinate, 
and can be expressed either in the r, u coordinates, or by using $ as a radial coor- 
dinate. 

Einstein's equations can be rewritten in terms of X{^) as: 

$' = -IvToXe-tZ-oc^'i*, (2.19) 
3 

yi' = (2.20) 
9 

Here, Wo > is the asymptotic value of the superpotential as $ — > — oo. It is related 
to the asymptotic AdS^ scale i by 

Wo = ^. (2.21) 
From these equations, the superpotential is related to X as: 

We deduce that fixing the function X{\) (hence the /3-function) specifies the super- 
potential (up to an overall multiplicative constant). Then, the equations of motion 
( |2.19 ) and ( p.2C ) determine the geometry up to two integration constants and an 
overall length scale i. As we show in Section 7, only one of the two integration con- 
stants is physical, the second one being an artifact of reparametrization invariance. 



2.3 Infrared properties of asymptotically AdS;^ backgrounds 

In the holographic approach to strongly coupled gauge theories, confinement at low 
energies is typically related to the termination of space-time at a certain point in 
the radial coordinate. In five- dimensional holography, with asymptotic AdS^ in the 
UV, this often implies the presence of a singularity in the bulk. We show here that, 
unless the IR is conformal, a curvature singularity is always present when we restrict 
ourselves to a two-derivative effective action. Specifically, we prove the following 
statement: 

Proposition: Consider any solution of 1 ^2.8^ ) such that exp A{r) i/r as r 
(with r > 0). Then, 

• The scale factor e^^"^^ is monotonically decreasing 

• There are only three possible, mutually exclusive IR behaviors: 

1. there is another asymptotic AdS^ region, at r ^ oo, where expy4(r) ~ 
i'/r, and i' < i (equality holds if and only if the space is exactly AdS^ 
everywhere); 
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2. there is a curvature singularity at some finite value of the radial coordinate, 
r = ro; 

3. there is a curvature singularity at r ^ oo, where the scale factor vanishes 
and the space-time shrinks to zero size. 

That the scale factor must be monotonic in this context is well known, and it is 
most clear in the u coordinates: the first equation in (|2.8|) implies that A"{u) < 0, 
therefore A'{u) must be monotonically decreasing. In the UV, A{u) ~ —u/i so for 
any u we must have 

A\u) < -l/£ < 0, Vm. (2.23) 

As a consequence, A{u) itself must be monotonically decreasing from the UV to the 
IR. This is a version of the holographic c-theorem [^. 

We now investigate possible IR behaviors. In conformal coordinates, the bound 
( p.23|) translates to : 



— e-^M > -. (2.24) 
dr i 

Suppose that the r coordinate extends to +oo. Then, either the l.h.s. of ( p.24| ) 

asymptotes to constant > £~^, or it asymptotes to infinity. In the former case, 
we obtain: 

f 

- ^ oo. (2.25) 



P' 
r ' 



This implies that the space-time is asymptotically AdS^ in the IR, with a smaller 
AdS^ radius i' . The gauge theory flows to an IR conformal fixed point, and is 
therefore not confining.^ 

If instead ^ e~'^'^'^^ — * +oo as r ^ +oo, then the curvature scalar diverges, as 
can be seen from its expression: 

R[r) = -e-'^ (l2A^ + Si) . (2.26) 

In this case, e"^"^*^'"-' diverges faster than r^, and A^ and A do not vanish faster than 
r~^, forcing eq. ( |2.26| ) to diverge as r — > oo. Moreover, the scale factor e"^'-''-' vanishes 
for large r, as claimed. 

There is another possibility, i.e. that the space-time ends at a finite value tq. 
This can happen because the scale factor shrinks to zero, or some of its derivatives 
diverge^ at tq. In either case eq. ( p^.26D indicates that we are in the presence of a 
curvature singularity at tq. 

These considerations were derived in the context of 5D Einstein-Dilaton gravity, 
but they are more general, because they follow only from the condition A"{u) < 0. 



^Five-dimensional holographic duals of the Bank-Zaks fixed points are in this class, pof . 

^We are always assuming that the space-time terminates due to some non-trivial dynamics, 
rather than because of some boundary at some otherwise regular point r = vq. This is in contrast 
with the original AdS/QCD approach which advocates an AdS space with an IR boundary. 
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This can be shown to be equivalent to the Null Energy Condition (NEC) (see e.g. 
[pp. Therefore the arguments of this subsection can be extended to any bulk field 
content, provided its stress tensor satisfies the NEC. 



3. Confining backgrounds 

Here we would like to characterize the backgrounds that exhibit confinement. By 
"confinement" we understand an area law behavior for the Wilson loop. Our anal- 
ysis allows a simple classification of confining background in terms of the metric, 
superpotential, or /3-function IR asymptotics. 

3.1 The Wilson loop test 

In this subsection we review the holographic computation of the Wilson Loop, 
12]. The potential energy E{L) of an external quark-antiquark pair separated by a 
distance L and evolved in time T, can be computed holographically by the action of 
a classical string embedded in the 5D space-time, with a single boundary which is a 
rectangular loop with sides L and T on the AdS^ boundary. We have, 

TE(L) = ^^G[Xl„(a,r)]. (3.1) 

Here Sng is the Nambu-Goto action evaluated on the world-sheet embedding X^^^{a, r) 
with minimum area: 

Sng = Tf J drda^/ - det gs, {gs)a/3 = igs)f.udaX^d^X'' , a,P = l,2 (3.2) 

with Tf = the fundamental string tension and {gs)fiu the bulk string frame 
metric. For a generic 5D metric of the form: 

{gs)iiudx^'dx'' = gss{s)ds'^ - gQo{s)dt^ + g\\ {s)dx^, (3.3) 

1^ showed that, for differentiable world-sheets, one can write: 

E{L) = Tff{sp)L - 2Tf ds^^^P{s)-P{sj.) 



(3.4) 



where the functions f{s),g{s) are defined as: 

f{s)=goois)g\i{s), g'^i-s) = goo{s)gssis) (3.5) 

and Si? is the turning point of the world-sheet in the bulk. Implicitly, sp depends on 
L through the relation: 

L = 2 r ds ^}^""] ^ ^ , (3.6) 
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where sb denotes the boundary. For large L, the second term in eq. ( |3.4|) is sub- 
leading. 

Expressions ( p.4|) and (|3.6| ) drastically simplify if we use conformal coordinates, 
s = r, 

(^5U(r) = e2^«(%,, As{r) = Air) + hir), f{r) = g{r) = e'^^^'\ (3.7) 



to obtain: 



L = 2 I dr , (3.^ 



/o ^g4As(r)-4As(rj.) _ ]_ 

In the neighborhood of r = the integral is finite, because the integrand behaves as 
g-2As(r) ^ g^j^j^ ^ ^3 small L. Around vp we may expand the denominator 
as: 

(3.9) 



^e4As(r)-4As(rf ) _ 1 ^4^'^ (r^r) (r^ - r) + 8A'^(rF)(rF - r)2 + . . . ' 

The integral is finite for generic rp and grows indefinitely as approaches a sta- 
tionary point of yls'(r), where A^ir^ = 0. This must correspond to a minimum 
since As starts decreasing at r = 0. In other words, if there exists such a stationary 
point r^, then 

rp ^ as L ^ oo. (3.10) 



The large L behavior of the quark-antiquark potential energy is thus (from ( |3.4| ) 

E{L) ~ Tfe^^^^'*^L (3.11) 

which exhibits an area law if and only if is finite at the minimum. From eq. 
( p.ll|) we read off the confining string tension as. 



T, = TfC^^^^'*^ (3.12) 

Notice that the finiteness of the string tension is not directly related to the value 
of the metric at the end of space, as sometimes assumed. Even if the space-time 
shrinks to zero-size at the singularity, which is the generic behavior of the Einstein's 
frame metric, this does not impede an area law: the string frame scale factor has 
a global minimum at a regular point in the bulk, and classical string world-sheets 
never probe the region of space beyond that point and never reach the singularity. 



Equation ( |3.11| ) captures the leading behavior of the quark-antiquark potential. 



In QCD the first subleading correction is the Liischer term, ~ 1/L. As shown in [14 



this term arises in some confining backgrounds (e.g. |T3[) from the first quantum 
corrections to the classical Wilson loop in It would be interesting to see if this 
is also the case in the models we are considering. 
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3.2 Confining IR asymptotics 

We are now ready to answer the question: which IR asymptotics give rise to confine- 
ment. 

Here we discuss a special class of metrics, that demonstrate particularly interest- 
ing features: namely the space-times with infinite range of the conformal coordinate, 
r G (0, oo). In Appendix ^ we give a complete discussion including other types of 
backgrounds. There, we also present the asymptotic values of some of the interesting 
quantities. The reader can find a summary of the classification in Table |l] at the end 
of this section. 

Consider a class of space-times whose Einstein frame metric has the form ( p.6| ), 
with the asymptotics: 

A(r) ^ -Cr" + . . . , r^oo, a, C > (3.13) 

up to generic subleading terms. Here, C~^^" = R is a length scale controlling the 
IR dynamics. 

The singularity is at r ^ oo, and the space-time shrinks to zero-size there®. To 
check whether the fundamental string is confining we need the string-frame scale 
factor, 

Asir) = A{r) + ^<l>(r). (3.14) 

As we have discussed in the previous subsection, confinement is equivalent to the 
existence of a minimum of the expression ( |3.14 ), where e^^ is non-zero. Due to the 
AdS UV asymptotics. As +oo as r — > 0. Therefore a necessary and sufficient 
condition for confinement is that As does not asymptote to — oo at the IR singularity, 
r +00''. 

The asymptotics of the dilaton can be obtained using the first equation in (|2.11|) : 

$(r) ~ -^A(r) + ^ log |i(r) | + $0. (3.15) 

Indeed, ( |3.15| ) solves eq. ( p.ll| ) up to a term proportional to {A/AY ~ ''^"^j regardless 
of the subleading behavior in ( pj.l3| ). 

Using ( p.l5|) we obtain the asymptotic form of the string frame scale factor (|3.14| ): 

As ~ -log|A(r)| r^^—^logr/R, ds^ ~ (-j {dr^ + r]i,dx'dx^) (3.16) 

Notice that the leading power-law term has canceled! Moreover the first surviving 
term is completely determined only by the leading power divergence of the Einstein 
frame scale factor. 



^The singularity is at a finite value ujr of the u coordinate. See appendix]^. 

''Since we are assuming that the singularity is at r — *■ cx), and 4> is monotonically increasing from 
(f> = —00 at r = 0, As cannot diverge to —00 at some finite r. Therefore, if there is a minimum for 
As, the string tension is certainly finite. 
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With the simple result ( |3.16 ), we can immediately determine which backgrounds 
lead to confinement: 



a > 1 =^ confinement: 

the string frame scale factor approaches +00 in the IR, thus it has a minimum 
at finite r. The special case a = 1 also leads to confinement. The minimum is 
reached as r — >■ 00, and the confining string tension is Tj lim^^oo exp[2y45'(r)].* 
Notice that when a = 1 the asymptotic geometry (in the string frame) is 5D 
Minkowski space-time with linear dilaton. 

a < 1 =^ no confinement: 

As asymptotes to —00 for large r, hence the confining string tension vanishes. 
It is easy to show that the same result applies if a = 0, and the scale factor 
A{r) goes to —00 slower than any power-law (e.g. logarithmically).^ 



We can relate the asymptotics ( p.l3 ) to the /9-function and to the superpotential, 
as follows: first we compute the X- variable, defined in eq ( p.l8|) , as a function of r, 



then, using eq. ( p.l5|) we can invert asymptotically the relation between $ and r and 
substitute it in the expression above. This gives: 



3a - 1 1 

1 + 7 T + • 

4 a log A 



A ^ 00. (3.17) 



We note that, generically, the point where A'g = corresponds to X = —1/2. In 
( |3.17|) , the point X = —1/2 is first reached at r*, and then at the singularity r = +00 



where A diverges. 



The asymptotic form of the superpotential is, from eq. ( |2.22| ): 



W{<^>) ~ $%re2*/3^ $ ^ +00. (3.18) 

Notice that in the leading asymptotics of the superpotential or of X(A) there is 
no trace of the dimensionfull constant C that controls the "steepness" of the warp 
factor in eq. ( |3.13 ). The appearance of the parameter R = C^^/" in the metric is 



the manifestation, in conformal coordinates, of the dynamical generation of the IR 
scale, as we will show explicitly in Section 7. It is fixed by the integration constants 
of Einstein's eqs, rather than by fundamental parameters appearing in W{X). 

The idea that some aspect of the geometry, which determines the IR scale, can 
be related to the integration constants rather than some a priori chosen parameter. 



®One could think of a situation where the string frame scale factor has multiple minima at r^, 
with nonzero values for exp[yl(ri)] (otherwise there would be a singularity at finite r.) In this case, 
the classical analysis implies that the string world-sheet has to stop at the minimum closest to the 
AdS boundary, and never knows about the existence of the others. However, quantum corrections 
may plausibly trigger the decay into other minima with lower tension. We don't consider this 
possibility any further, and we will always treat backgrounds with a single minimum of As- 



As discussed in Section 2.3 A{r) cannot asymptote to a finite constant. 
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was already present in the "braneless approach" to AdS /QCD of [15|. As we will 



discuss in section |6.6.1| however, the spectral properties of the background analyzed 
in |]I5[ suffer from some pathologies, that make it conceptually equivalent to models 
with a hard IR cutoff, in which some additional, arbitrary boundary conditions in 
the IR must be supplied. 

We can also relax the requirement that A{r) grows as a simple power-law, since 
from eq. (|3.16|) we see that all that is needed for confinement is the condition 
logl^l > asymptotically. This is true for any function A{r) whose asymptotics is 
bounded above and below as: 



Cir"^ < -A{r) < C2r"^ ai,2 > 1, Ci,2 > (3.19) 
3.3 General confinement criteria 

In Appendix ^ we analyze also the backgrounds where the singularity is at finite 
r = tq. They always exhibit area law. The analysis in the previous section, together 
with Appendix 0, allows us to formulate a general criterion for confinement in 5D 
holographic models: 

General criterion for confinement (geometric version) 

A geometry that shrinks to zero size in the IR is dual to a confining 4D theory if 
and only if the Einstein metric in conformal coordinates vanishes as ( or faster than ) 
e"*"*" as r ^ oo, for some C > 0. 

(It is understood here that a metric vanishing at finite r = Tq also satisfies the above 
condition.) 



Comparing the superpotentials found in all the examples studied in Appendix 
0, eqs. ( |A.21|) , ( |A.36|) , (|A.50|) , and ( |A.65|) , we see that one can treat simultaneously 
all cases by using the following parametrization for large A : 



W{\) ~ {\og\f'^\^, /5(A) = 3AX(A) ~ -^A 



P 1 
2" log A 



(3.20) 



where P and Q are real numbers. This implies for the Einstein and string frame 
dilaton potentials: 



V{^) ~ (log A)^A2«, Vs{^) ~ (log A)^A(2Q-4/3) 



(3.21) 



An equivalent characterization of the confining backgrounds is: 



General criterion for confinement (superpotential) 

A 5D background is dual to a confining theory if the superpotential grows as (or faster 
than) (log A)'^/^A^/^ as A oo for some P >0. 
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The relation between parameters P and a (appearing in (|3.13 ), is given in Table 
One can also relate the IR properties directly to the large A asymptotics of the 



/3-function. Computing X(A) = /?(A)/(3A) from the superpotential via eq. (|2.22| ), 
one obtains the following form of the same criterion: 

General criterion for confinement (/^-function) 

A 5D background is dual to a confining theory if and only if 

lim (X(A) + M logA = if, -oo<K<0 (3.22) 



A^oo \ 2 

In the above form^° the condition for confinement does not make any explicit refer- 
ence to any coordinate system. Yet, we can determine the geometry if we specify K. 
In particular: 

1. K = —oo: the scale factor vanishes at some finite ro, not faster than a power- 
law. 

2. — oo < K < —3/8: the scale factor vanishes at some finite ro faster than any 
power-law. 

3. —3/8 < i^' < 0: the scale factor vanishes as r cxd faster than e"'"''^^' for 
some e > 0. 

4. if = 0: the scale factor vanishes as r ^ oo as e"*"^ (or faster), but slower than 
g-Cri+" for any e > 0. 

The borderline case, K = 3/8, is certainly confining (by continuity), but whether 
or not the singularity is at finite r depends on the subleading terms. When K is 
finite, we can relate it to the parameters Q and P appearing in the superpotential: 
if if < oo, then Q = 2/3 and P = —8K/3. The classification of the various possible 
IR asymptotics in terms of their confining properties is summarized in Table |l|. 

We note that, if we classify the backgrounds in terms of P and Q, our analysis 
covers the entire range of these parameters. As a result, our classification is not 
limited to superpotentials that behave asymptotically as ( |3.20| ), but it also applies 
to any superpotential that for large A is bounded between two functions of the form 
( |3.20|) , for two appropriate pairs {Qi,Pi) and {Q2,P2)- 

For most of the confining backgrounds, although the space-time is singular in the 
Einstein frame, the string frame geometry is regular for large r (see Appendix ^). 
In fact, in these situations, all curvature invariants vanish for large r. The dilaton 



"'^^We are assuming that the hmit exists, and that the /3-function does not oscillate infinitely many 
times across — 3/(2A) as A — > oo. This possibility seems remote from a physical point of view. 
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r e (0, oo) 


r e (0, ro) 


A{r) 


—7 log r 




-C(ro - r)-" 


5 log(ro - r) 






< a < 1 


a > 1 






Confining 


No 


No 


Yes 


Yes 


Yes 














Q 




2 
3 


2 
3 


2 
3 




p 


arbitrary 


^ < 

a 


^ e [0, 1) 


a 


arbitrary 


K 


> 


> 


-i^e(-|,o] 


-|^e(-oo,|) 


—00 



Table 1: Summary of confining asymptotics. As required by the NEC, tlie parameters 
a, d, 7, 5, C are all assumed positive, and 7 > 1 



however diverges. Therefore, in the string frame the singularity manifests itself as a 
strong coupling region in a weakly curved space-time. 

Interestingly, as discussed in the previous subsection, string world-sheets do not 
probe the strong couphng region, at least classically. This is because the geodesic 
surfaces ending on the AdS boundary do not stretch beyond the minimum of the 
scale factor. At that point, the t'Hooft coupling may be of order one, but the string 
coupling Qs = \/Nc is still small. This can be intuitively attributed to the fact that 
the string tries to stay away from the region where the metric becomes large, since 
this would generate a larger world-sheet area. 

Therefore, singular confining backgrounds have generically the property that the 
singularity is repulsive, i.e. only highly excited states can probe it. This will also 
be reflected in the analysis of the particle spectrum, in the next subsection. This 
consideration makes our conclusions more robust, since they are insensitive to the 
region near the singularity, where quantum effects may become important. As the 
classical string world-sheet never probes the strong curvature region, a semiclassical 
analysis is reliable. 

One could also worry that a direct coupling of the dilaton to the world-sheet 
curvature scalar could spoil this analysis. This is not so, as shown in Appendix C. 

3.4 Magnetic cheirge screening 

In confining theories, one expects the dual magnetic gauge group to be Higgsed, lead- 
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ing to a screening of the magnetic charges. In our setup, magnetic monopoles can be 
described as the endpoints of Dl-branes. Therefore the calculation of the monopole- 
antimonopole potential proceeds exactly like the one for the quark- ant iquark poten- 
tial, with a D-string replacing the fundamental string. In this section we discuss the 
case of infinite range backgrounds, leaving the finite case to Appendix |B[ 
The D-string action is^^ 



Sd = Td^J d^ie '^^J-det g^p, (3.23) 

where ga/3 is the induced metric on the world-sheet and the target space metric is in 
the string frame. We work in the conformal coordinates, 

ds^ = e^^^ {dr^ + r]ijdx'dx^) , (3.24) 

and reabsorb the factor of the dilaton in the conformal factor of the target space 
metric, reducing the problem to a string with Nambu-Goto action propagating in a 
target space with an effective metric: 

ds^ = e^^^^'-) {dr^ + riijdx'dx^) , Ad = As-^ = A + ^. (3.25) 

The properties of the string embedding can then be deduced using the same tech- 
niques as in the previous subsections. 

For large L the energy of this a configuration is simply given by 

E = e2^^('^^(^))L + ■ • • , (3.26) 

where rp{L) is the bulk position of turning point of the worldsheed with length L 
on the boundary. The relation between rp and L is given by a formula similar to eq. 
(|3.8|), with the substitution As — > Ap,. 

To avoid the magnetic charge confinement, it must be that the scale factor 
of the "D-string frame" metric, eq. ( ^.25| ), vanishes at the IR singularity. 

In the confining backgrounds of section |3.2| , with large r asymptotics ( p.l3| )- 
( |3.15D and a > 1, the scale factor e^-° in eq. ( ^.25| ) does indeed vanish as r — >• cxo; 



the magnetic string tension is zero and the magnetic charges are not confined. The 
question remains, whether they feel an inverse power-law potential or they are truly 
screened in which case the potential falls-off exponentially or faster. Below, we show 
that the latter holds for the backgrounds under consideration. 



"'^^One might expect that the magnetic and electric quark potentials are related by A 1/A 
duality, in the Einstein frame. This is the case for gauge theories that are deformations of A/" = 4 
sYM and others that descend from ten dimensions but not in general. In the Einstein frame the 
electric string action is proportional to A 0-2 while the magnetic one to X°-^ . For D = 10, these 
factors are inverses of each other, but not for D — 5 relevant here. 
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In order to answer this question, one has to study the potential energy (|3.26|) for 
large L: to do this, one has to invert asymptotically the relation between rp and L 
from the D-string analog of eq. ( |3.8| ), and insert it into eq. ( |3.26| ). 

The asymptotic form of the D-string metric is 



a > 1, 



where we are restrict to the confining case. We must evaluate 



rF) 



rp 



dr 



- 1 



1/2- 



(3.27) 



(3.28) 



By assumption, there are no other singularities of exp[y4£)] for any finite r, and we 
assume that there are no other local extrema. Thus, the only region in which L{rp) 
could diverge is r^r — oo. 

We show below that, for a > 1, L^^^rp) is finite in this limit. We first divide 
the integration range in two regions, < r < ri, ri < r < rp, such that in the second 
region the asymptotic form of the scale factor ( 3.27|) holds. Consider the integral in 
the first region: 



-ri rri g4(Ai3{rF)-Ai3{ri)) 

dr- 



Jo [e^(AD{rF)-Ao{r)) _ l^^/^ Jq [e4(AB(r)-Ao(ri)) _ g4(Ao(rf )-Az,(ri))] 1/2 

< / 777 = Lin). (3.29) 

The inequality follows from our assumption that Ad is monotonically decreasing. 
Since L{ri) is finite for finite ri, and explAA^irp)] ^ as rp ^ oo, the r.h.s 
vanishes in this limit. Therefore, for large rp the dominant contribution to L{rp) 
comes from the region r > ri. 

To analyze the behavior of the integral over the asymptotic region, consider first 
the case a = 1. We have: 

and 

1 /"^°° dv 

lim L^^\rp) = — = Lmax < +oo. (3.31) 
Next consider a > 1. For large r < rp, the following inequality holds: 

r'^-r'' > r'^-\rp-r). (3.32) 
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It follows that 



1/2 



(3.33) 

which implies that for a > 1 

lim L^'^\rF) = 0. (3.34) 

rp^+oo 

We showed that L^°'\rF) cannot be larger than a maximum value Lmlx, which is 
reached at +oo if a = 1, and at some finite Tmax if a > 1. Therefore two monopoles 
at a distance larger than L^ax cannot be connected by a smooth world-sheet. In 
this case, the configuration that minimizes the action consists of two straight lines 
separated by a distance L and joined by a line at constant r = oo. This configuration 
has the same energy as the one with two straight lines only as the contribution from 
the piece at the singularity vanishes^^. Therefore for L > L^ax the monopoles are 
non-interacting. This shows that in the backgrounds with an infinite range of r, the 
magnetic charges are screened. 

The finite tq case is discussed in Appendix 0, where we show that the monopole 
charges are screened, except in backgrounds with power-law decay expA ~ {tq — r)^ 
with 5 < 1/15. This case falls into the range < 5 < 1, which turns out to 
be problematic also for other reasons as we show in the discussion of the particle 
spectrum in Section 

3.4.1 Absence of screening in hard- wall models 

In the simplest models proposed as a holographic description of chiral dynamics of 
QCD 1^, 1^, the space-time ends at an IR boundary before any singularity. According 
to our discussion in this section, one finds linear confinement both for the electric and 
the magnetic charges. This is contrary to the expectations from the gauge theory 
dynamics. In fact, the computation of the magnetic string Wilson loop is exactly the 
same as that of the electric one, since the wall has the same effect on both objects. 



This was computed for the cut-off AdSs background in e.g. |]T6|, where the expected 
area law was found. 



4. The QCD vacuum energy 

An interesting question in YM theory concerns the value of the vacuum energy, and 
this is closely related to the so called gluon condensate, (Tr[F^]). Typically they 

^^One should take this argument with a grain of salt. This is because, unlike the configuration 
that stretches only up to rp, this configuration falls into the singularity, hence one should worry 
about various string and quantum corrections to the classical solution. At any rate, our final 
statement about the magnetic screening is valid as existence of an Lmax is sufficient for that. 
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are UV divergent. Of course one would try to renormalize them by subtracting the 
divergences. We do not know of an unambiguous way to define them beyond per- 
turbation theory^^. Indeed, once the divergences are subtracted one might as well 
subtract also the finite piece. We also stress that the gluon condensate is also de- 
fined in a semi-phenomenological fashion via the SVZ sum rules. Therefore, without 
quarks, it is not obvious how to define it^"^. 

Because of this we will calculate the divergent (holographic) full vacuum energy 
to leading order in l/Nf.. To do this, we will introduce the usual UV cutoff near the 
AdSs boundary and will compute the Euclidean action of the vacuum QCD solution. 
The Gibbons-Hawking boundary term is important in this calculation. 



— Se + Sgh 



Sgh = 2M^Nl 



R - -{d^f + V{^) 



d^xVh K 



dM 



Evaluated on a solution the Einstein action is 



Se = ^M^N^ J d^x^ V{^) = ^M^N^Vi ° dr e^^ V{^) 

(e^^i) = 2M^N^V4 



(4.1) 
(4.2) 

(4.3) 
(4.4) 



dr 



where V4 is the space-time volume and we introduced a IR cutoff e in the bulk. For 
all our confining backgrounds, the contribution from the singularity vanishes auto- 
matically. This is a good consistency check of the procedure as only the contribution 
from the UV boundary should survive. We therefore obtain 

Se = -2M'^N^V4 e^^^^Mje) (4.5) 

For the GH term, the trace of the extrinsic curvature is = Ae~^A and therefore 



(4.6) 



Putting everything together we obtain for the vacuum energy density (Euclidean 
action divided by the space-time volume) 



(4.7) 



^^See however for another discussion. 

-'^^This is unhke the CP-odd condensate {Tr[F A F]) that we wiU calculate in section ^ The 
reason is that the operator Tr[F A F] does not need a holographic renormalization. 
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Note that this is negative as A < 0. We can re-express the result in terms of the 
cutoff energy scale, Auv 

A{e) = log Auv 

Then the bare vacuum energy density satisfies 

91ogAc/y 9 9A2 9 ^ ^ ^ ' 

Another (related) equation determines the coupling dependence 

aiog^QCD 4 (3{X) 



dX{Auv) /3(A) 4A2 
5. The axion background 



(4.9) 



The axion a is dual to the instanton density Tr[FAF]. In particular its UV boundary 
value is the UV value of the QCD ^-angle. Moreover, its profile a(r) in the vacuum 
solution may be interpreted as the "running" 6'-angle in analogy with the dilaton, 
that we interpret as the running coupling constant . Such an interpretation should 
be qualified, as it may suggest the false impression that UV divergences renormalize 
the 6 angle in QCD. We will return to this later. 

The question of the 9 dependence of large Nc QCD and the associated rj' prob- 
lem has led to several advances that culminated with the Witten-Veneziano solution, 
[PH , p^ . It states that although naively the 9 dependence is expected to be non- 
perturbative, at large this is not so. It enters at order in YM theory. It 

generates a 6'-depended vacuum energy that scales as and provides the correct 
mass (of order 1/Nc) to the rj'. Such expectations have been verified in the holo- 
graphic realization of a four-dimensional confining gauge theory based on D4 branes. 

Here we analyze the structure of the background solution for the axion in five 
dimensions relevant for pure YM theory. The action in the Einstein frame and the 
corresponding equation of motion are: 

= ^ y d'xy^ Z{X) {d,af , -^d, [Z{X)^g'^''d,] a = (5.1) 

where Z{X) captures a part of the a' corrections. It was shown in appendix B.l of 
[0] that Z{X) depends in general on the 't Hooft coupling A. In perturbative string 
theory and to leading order in a', Z{X) = X^. However, as was the case for the 
potential, we expect a constant leading term also here, 

Z{X) = Za + 0{X') , A^O (5.2) 

The axion field equation is to be solved on a given metric and dilaton background, 
i.e. we neglect the backreaction of the axion 
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For a radially dependent axion the equation becomes 

'd+ (3A + {dxlogZ)\^ d = (5.3) 
This equation can be integrated once as 

C e"^^ 

The equation (|5.3| ) has two independent solutions. One is a constant, foir) = 
constant. The other /i(r) can be obtained by integrating ( p.4| ) and choosing the 
initial conditions so that it vanishes at the boundary r = 0: 

where we divided by the AdS length so that the function is dimensionless. A first 
observation is that the function /i(r) is strictly increasing. 

Since near the boundary, Z = + ■ ■ ■ , e* = — ^Qg^.^^) + ■ ■ ■ and = ^ + • • • 
we obtain 

r „ / 1 

lim /i(r 



1 + 



log(rA) 



(5.6) 



where we chose an arbitrary normalization for this solution. This solution is the one 
normalizable in the UV. 

The constant solution should be related to the UV value of the ^-angle as 

/o = euv + 2TTk , keZ (5.7) 

The different values of the integer k correspond to an infinite number of vacua, [p^] . 
Such vacua exhibit in general oblique confinement and they are degenerate to leading 
order in the 1/Nc expansion. As it is known and we will also see it explicitly below, 
their degeneracy is lifted to the next order. They can be separated by domain walls 
that are the D2 branes. 

The full background solution therefore reads 

a(r) =^t/y + 27rfc + C /i(r) (5.8) 

where we take by convention 6uv ^ [0,27r). The coefficient C is proportional to the 
expectation value of the QCD instanton density operator in the QCD vacuum. Using 
the precise holographic formula and ( |5.6| ) we obtain 

AZJ' {FAF) _ {FAF) 
^ " (2A - 4) 327r2 " ^'^^ 327r2 ^^■^> 

Substituting the solution in the effective action we obtain, for the (^-dependent) 
energy per unit three-volume, the following boundary terms 

^{Ojjv) = — dr^Zmdaf = —e'^Z{\)aa = —C a(r) (5.10) 
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where we have used the equations of motion to write the on-shell action as a boundary 
term, tq is the position of the singularity in the IR. It may be finite or infinite, as 
discussed in the previous sections. 

We expect that the only contribution to the 9 dependent vacuum energy should 
come from the UV boundary. The reason is that there should be only one boundary 
in the theory. The presence of a second boundary would imply that the holographic 
dynamics of the theory is incomplete. Therefore, we should not expect a contribution 
from r = tq. In order for this to be true, the axion should vanish at the singularity.^^ 
We must therefore have, 

Eieuv) = ^C{9uv + 27rk) , air^) = 9uv + 27rk + C Mro) = (5.11) 
Solving the IR equation assuming /i(ro) 7^ we obtain 



2i /i(ro) ' euv + 27rk 



Mr) 



ro dr 



(5.12) 

where the minimum on k (that we denote by /cq) is obtained in order to choose the 
k-vacuum that minimizes the energy. /^From ( ^.12| ) we can extract the topological 
susceptibility as 

Jo e3'4z(A) 

We have obtained the expected quadratic and non-analytic behavior for E{9). 
We have also determined that the instanton vacuum condensate is non-zero 

{FAF) _ Ouv + 27rfco _ Ouy + 27rA;o 
32vr2 Z.£Vi(ro) ^ ^ ' ^ 

Notice, that the constant Z(^) drops out of all quantities of interest associated with 
the axion, except the topological susceptibihty. Moreover, the condensate is finite 
without a renormalization of the FAF operator. This is in accordance with lattice 
results . 



We also observe a very interesting corollary: // we view a{r) as an effective 
9 -angle (in analogy with the t' Hooft coupling) it vanishes in the IR! . 

We now study the dimensionless constant /i(ro) that is proportional to the in- 
verse of the topological vacuum susceptibility 



"'^^ Allowing any other value at the singularity, a(ro) = we obtain for the vacuum energy density, 
£ ^ (dyy — ao)^ in contradiction with all large- A^c expectations, [^, p2[ . 
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The integrand is a positive function as Z{X) is multiplying the axion kinetic energy 
and is therefore expected to be non-negative. Moreover we do not expect the in- 
tegrand to diverge at a point before the singularity rg, as vanishes only at tq, 
and Z{X) is also not expected to vanish. Therefore, the only potential pathological 
behavior is a divergence at tq. 

To study the region around the singularity we will have to study the two cases 
(ro finite or infinite) separately. 

• We first consider the IR asymptotics in the infinite range case, namely the 
singularity at r = oo. From section |3.2|, for large r and in the Einstein frame: 



logA = -Cr'^ + --- 
We also assume that for large A, 

Z{X) ^ + 

Then: 

(1) i{dy^2 



A = -Cr"" + ■ 



oo 



(5.16) 
(5.17) 



/i(ro 



oo 



J - exp[— (2-rf)r- + 



(5.1^ 



In order for this not to diverge, we ask d > 2. In this case the low energy 
asymptotics of the axion are 



a r 



1 



9uv + 27rA;o 
1 



a/i(oo) 



a/i(oo 
3{d-2)C 



/i(oo 
3{d-2)C 



drexp —-{d — 2)Cr'^ 



-,|(d-2)Cr^ 
a I 



(5.19) 



--1 



r" ^ exp 



--[d-2)Cr- 



^K-^-s) (logE)"^ 



where in the last expression we have replaced the radial variable with the 
holographic energy using ( |2.15|) . 

(2) For d = 2, 



fiiro 



oo 



dr r-i("-i) -I- 



(5.20) 



In order to obtain a finite result, a > 5/3. This is stronger than the confinement 
condition a > 1. The low energy asymptotics of the axion are 



a(r) 



Ouv + '^T^ko /i(oo) 



dr r 2 



26^ 



uv 



(3a-5)/i(oo) 



(3a-5) 
f 2 ^ 



(3a-5) 

-logE) 
(5.21) 



and the effective 6'-angle grows slowly in the IR. However, as it is shown in 
section |6.5.3| , in order for the O"*"" glueballs to have a discrete spectrum, we 
must demand d > 2 and therefore this case seems not relevant for QCD. 
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Similar remarks apply to confining backgrounds with tq finite. In particular 
/i(ro) is finite ii d >2. When d > 2 then at low energy 



e{E) ~ (5.22) 

while for d = 2, the low energy running is by the inverse power of the logarithm 
of the energy. 

A comment is in order here concerning the relevance of higher-derivative cor- 
rections to the IR asymptotics of the solution. Such corrections can be obtained by 
substituting Z{X) Z{R, (9A)^, A). This is indeed the most general form as higher 
powers of {daY are suppressed by extra powers of l/Nc. All arguments made above 
go through with this more general kinetic function. 

The oblique confinement vacua, labeled by k, are long-lived in the laxge-Nc limit, 
with lifetimes that scale as 0{Nc). The D2 brane, discussed in [|^ is a domain 
wall separating two such consecutive vacua, as k jumps by ±1 when crossing the D2 
domain walls. 

5.1 Screening of CP violation in the IR 

The essence of the strong CP problem lies in the fact that a non-zero ^-parameter 
in QCD breaks CP (except at 6' = tt) and provides a non-trivial contribution to the 
neutron dipole moment. The stringent experimental limits on this quantity constrain 
6 to be very small, {-< 10^^) [Q. This is known as the strong CP problem: why is 
6 so small in nature? 

In pure YM, the case that we are studying here, the nature of the strong CP 
problem changes, as there are no quarks, and no baryons. However, the issue is how 
strong are the CP violating effects of a nontrivial 6uv in observable data. For a pure 
gauge theory, a CP-violating effect can be the decay of an (excited) 0"*"" glueball to 
0"*"^ glueballs. How important such effects are at low energy depends on the way 
the axion solution behaves at all radial distances. The behavior of the radial axion 
profile is shown in figures |l| and |[ 

It is tempting to think of the radial axion solution as a "running 6'-angle" in 
analogy with the similar intuition concerning the dilaton that has been justified 
quantitatively in numerous holographic setups. Here however, such an interpretation 
needs to be qualified, as there are strong indications from lattice P3, that 6 as 



a coupling in the bare Lagrangian does not receive UV-singular corrections, that 
would force it to renormalize. That being said, it is direct on the other hand in the 
holographic context to see that parity violations at low energy although proportional 
to 6uv, have numerical coefficients that are due to this radial change of the axion 
and such coefficients can be small. 

It is well known that in theories where non-perturbative corrections can be con- 
trolled, that 6eff at low energy receives finite corrections and is different from Buy 
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A simple example of this are M = 2 gauge theories in four-dimensions. There, the 



effective field theory in the Coulomb branch can be exactly solved, [^], and exact 
effective superpotential calculated. Its real part is the effective 6'-angle, which does 
receive instanton corrections. Although here the effective theory is abelian the ef- 
fective ^-angle can affect low energy physics as some monopoles and dyons can be 
light. In large- iVc YM, instantons are apparently suppressed exponentially, but as 
was argued in early works, ||2T|, such non-perturbative effects can affect 6'-related 
physics at leading orders in We therefore advocate that the radial change of 

the axion can be interpreted as a finite renormalization of the 6'-angle in QCD. Of 
course the quantitatively precise statement is that the axion solution must be used 
in accordance with the standard rules of the gauge theory/gravity correspondence to 
calculate physical quantities. 
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Figure 1: An example of the axion profile (normalized to one in the UV) as a function of 
energy, in one of the explicit cases we treat numerically in Section ^. The energy scale is 
in MeV, and it is normalized to match the mass of the lowest scalar glueball from lattice 
data, mo = 1475MeV. The axion kinetic function is taken as Z{X) = Za{l + CaA'^), with 
Ca = 100 while it does not depend on the value of Za- The vertical dashed line corresponds 
to Aqcd as defined in eq. ( 7.10| ). In this particular case A = 290MeV. 



The discussion above as well as figures and ^ suggest that the CP-violating 
effects of 9uv are screened at sufficiently low energy, whatever the UV value of the 
^-parameter. Moreover, the IR vanishing of the axion is power-like as we have shown 
above. As we argue in section |6.5.3| , the expected value of the exponent d in ( ^.17| ) 



is d = 4 from parity independence of asymptotic glueball spectra. 

There have been attempts to calculate QCD data at small values of Ouv^ on the 
lattice. Results in particular exist on the topological susceptibility as well as the 
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Figure 2: A detail showing the different axion profiles for different values of Ca- The 
values are Ca = 0.1 (dashed line), Ca = 10 (dotted line) and Ca = 100 (solid line). 



^-dependence of the glueball spectrum, Such calculations can in principle be 

confronted with similar calculations using the present framework, but we will not do 
this here. 

The relevance of this discussion in the case where quarks are present is more 
complex and we will not attempt it here although we will show several of the relevant 
ramifications in section |6!7. 



6. The particle spectrum 

In gauge-gravity dualities, the particle spectrum of the 4D theory is obtained from 
the spectrum of fluctuations of the bulk fields around the background. Reviews of 
glueball spectra calculations in the holographic context can be found in ||2^. In 



this section we first give a general overview of the spectra of various particle species 
(glueballs and mesons). Then, in Section |], we compute numerically the glueball 
spectrum of some concrete backgrounds that exhibit asymptotic freedom in the UV 
and confinement in the IR. The main results of this analysis can be summarized as 
follows: 

1. In the previous section we showed that in order for the theory to confine, the 
Einstein frame scale factor must vanish at least as fast as exp[— Cr°] with 
a > 1, C > 0. Remarkably, this is the same condition one obtains from the 
requirement of mass-gap in the spectrum. Therefore, using holography, we can 
directly relate the existence of a confining string with the existence of a mass 
gap. 
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2. Among the class of confining backgrounds we liave considered, we find examples 
that exhibit an asymptotic "linear" mass spectrum, ~ n. 

In this section we are mostly interested in confining backgrounds in which the 
scale factor exhibits exponential decay at r ^ oo; in the last subsection we briefly 
discuss the backgrounds with a singularity at finite r. 

6.1 General properties of the spectra 

Here we discuss the spectrum from a general point of view and leave the details and 
comparison with the lattice results to Section p. We mostly work in the conformal 
frame, where the properties of the spectrum are more transparent. The spectrum of 
particles up to spin 2 is determined by the fluctuation equations of the various bulk 
fields in the solution. Typically, one can identify fluctuations ^(r, x*) with a diagonal 
kinetic term and a quadratic action of the form 

S[i\ ~ j drd'x e^^M + m? + M\r)e] , (6.1) 

where B{r) and M^(r) are functions depending on the background and on the type 
of fluctuation in question. 

The linearized field equation reads: 

^ + 2Bi + - M\r)^ = 0. (6.2) 

To look for 4D mass eigenstates, the standard procedure is to write: 

e(r,x) = ar)e'\x), n^^'){x) = m'^('\x). (6.3) 

Then, eq. ( |6.2D can be put into a Schrodinger form by defining a wave-function 
associated to the fluctuation ^, 



Eq. (p.2|) becomes 



^(r) = e-^(''V(r) . (6.4) 



dr"^ 



if) + V{r)ip = m ijj, (6.5) 



with the potential given by. 



, d^B f dBV , , , 

V{t) = --+( — ] +M\r). (6.6) 



dr'^ \ dr 

The Schrodinger equation ( |6.5| ) is to be solved on the space of square-integrable 
functions ip{r), as can be seen inserting (|6.4|) into the quadratic action: the kinetic 
term of a given 4D mode ^^^\x) reads: 



rfre2^W|e(r)|2 ) / d^x {d^^^^\x)) = ( / dr\ij{r)\^ ) / d'x {d^^^^\x)) . 

(6.7) 
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Requiring finiteness of the kinetic term leads to 



drlijj {r)\^ < oo. (6.8) 

Therefore, in these coordinates, the problem of finding the spectrum translates 
into a standard quantum mechanical problem. The general features of the spec- 
trum can be inferred from the properties of the effective Schrodinger potential (|6.6| ). 
Given the functions B{r) and M(r) we can obtain useful information without finding 
explicit solutions. 

In the case we are mostly interested in, i.e. the infinite-range case, a number of 
interesting properties of the spectrum can be derived in full generality. 

6.2 Existence of a mass gap 

Consider first the effective potential in the asymptotically AdS^ region, r ~ 0. There, 
the potential behaves universally, since B{r) ~ 3/2 A (r) in the UV for all kinds of 
fluctuations: 



i^^+oo. (6.9) 



15 1 
4 r 

Next, notice that the equation ( |6.5|) can be written as: 

{P^P + M^{r))tlj = m'^ij, P={-dr + B{r)) (6.10) 
Taking into account also the behavior near r = 0, the following statements hold: 

1. if M'^ijr) > the spectrum is non-negative 

2. If the space-time extends to r = -\-oo, and if moreover V{r) does not vanish as 
r oo, then there is a mass gap. 

The second statement follows from the fact that, if the singularity is at r ^ +oo, 
there can be no normalizable zero- mode solutions to eq. ( |6.5|) with M^(r) = 0. 
Indeed, the m = solutions can be found exactly: 

(r) = e^('-), = e^('-) dr'g-^^^^^'') (6.11) 

Jo 

It is easy to see that ifj^^^ is not normalizable in the UV, since there B ~ —3/2 logr. 
On the other hand, ipQ is not normalizable in the IR, no matter what is the behavior 



of i?(r) as r 00. See |TT| for a detailed discussion. Since there are no normalizable 
zero-modes, the only obstacle to having a mass gap would be a continuous spectrum 
starting at zero. This cannot be the case if V{r) does not vanish as r — 00. 

If the singularity is at a finite r = tq, this argument still holds if the scale factor 
vanishes at least as (ro — rY, with 6 > 1. Then, again one can show that there are 
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no zero-modes. For < 5 < 1 things are more subtle, and we will leave this case for 
separate discussion in Section |6.6.1 . 

For the various particle types we analyze (vector mesons, and glueballs of spins 
up to 2) we will see that property 1 always holds, in fact = 0. Moreover, for all 
particles we consider, it turns out that the function B{r) has the same IR asymptotics 
as A{r). In the backgrounds with infinite r range, as r ^ oo: 

^(')--{^y^ (6.12) 

therefore 

V{r) = B^{r) + B{r) ~ R'"^ (^)^^'' • ^^•^2) 

We see that the mass gap condition is a > 1. This is the same condition we found 
independently for quark confinement. If we require a > 1 strictly, we moreover obtain 
a purely discrete spectrum, since then V{r) —>■ +oo for large r. If a = 1 the spectrum 
becomes continuous for > V{r — > oo). 



6.3 Large n mass asymptotics and linear confinement 



In the confining backgrounds, where the potential behaves as in eq. (|6.13| ) for large r 
and as ( |6.9| ) for small r, the large eigenvalue asymptotics of eq. ( |6.5| ) may be obtained 
through the WKB approximation: the quantization condition is approximately given 
by the quantization of the action integral: 



mr 



r2 



(6.14) 



where r\ and ri are the turning points. For large m^, r\ ~ 0, and {ri/Rf'^"' 
R^m?, so we can write: 



UTT = m 



i?(m,i?)i/("-l) 



(6.15) 



Assuming ^ V{r) in the intermediate region, the second term under the square 
root becomes relevant only when V{r) takes its asymptotic form. We can therefore 
write 



W 1 — 

where A = R~^. For large n: 



"1 1 



R/ mR 



dr 



A 

m ~ A n Q . 



(6.16) 
(6.17) 
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In particular we have "linear confinement" (m^ ~ ra) if a = 2}^. For a — »• oo 
the spectrum looks similar to the one of a "particle in a box" potential, ~ n^, 
characteristic of "hard wall" models and more generically of any background with 



finite Tq (see Section p.6|) . 



6.4 Universal asymptotic mass ratios 

Here we derive some general properties of the glueball spectrum that are independent 
of the specific potential chosen. In this section, we consider the backgrounds where 
X —1/2 at the singularity. As we have seen, this is generic in confining back- 
grounds with singularity at r = +oo. The function B in (|6.2|) generally asymptotes 
to 

B ^Bf log{X), as A^oo, (6.18) 

where the coefficient Bf depends on the type of particle. It is essentially determined 
by the normalization of the kinetic term in the effective action of the specific type of 
background fluctuation that correspond to the particle in question. 

One can also write down the effective Schrodinger potential (|6^ ) using A as a 
coordinate, 

. ^..ir ^(*-^.) . (3^ . 1 - f . . (Af ) 

(6.19) 



From ( |6.19D we observe that Vg in the IR asymptotes to. 



K - ^e'^^'^'^r'^(^~''')B' > 0. (6.20) 

The exponential depends on the specified /3-function of the gauge theory. However 
the constant Bf is universal for a given type of particle, i.e. it is independent of 
the specified running of the gauge coupling. For example Bf is 1 both for the 0^"*" 
glueballs and the 2++ glueballs, i.e. 

K(0++)(r) ^ 



K(2++)(r) 



r^+00. (6.21) 



"'^^A dilaton and/or a warp factor A{r) behaving as for large r, were advocated in |29[ , in order 
to obtain a linear spectrum for mesons. In that work, the authors suggest an AdS^ space-time 
together with a dilaton with asymptotics. Such a background is sometimes referred to as a "soft 
wall" model, and has been used to compute meson-related quantities (see e.g. for recent work). 
A similar approach was adopted in to treat both baryons and mesons. In that work it is the 
scale factor, rather than the dilaton, that drives the IR asymptotics. 

We should stress that those backgrounds, unlike the ones we study here, are not obtained as 
solutions of any set of field equations. For example, from our previous discussion it is apparent 
that, if the dilaton grows as in the IR, its backreaction is such that the space-time cannot be close 
to AdS5 for large r, independently of the form of the dilaton potential. Moreover, as we discuss in 
Section 5.4, the dynamics of mesons could be described by a different mechanism |^ which does 
not necessarily require a — 2 for a linear meson spectrum. 
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This means that the glueballs have a spectrum whose slope is independent of 
their spin for large mass: 

This fits nicely with the semi-classical string models (see e.g. [Q) for the glueballs 
that predict 

2 

fry 

= 2n + J + c, (6.23) 

where aa is the adjoint string tension, J is the angular momentum and c is some 
number of order 1. Our finding ( |6.22| ) is in accord with the general prediction of such 



models that the adjoint string tension is universal for glueballs with different spin. 
Next, we move to specific analysis of the spectra of different species of glueballs. 

6.5 Glueball spectra 

At the lowest mass level the bulk theory contains the dilaton $, the metric g^^, 
and the axion a. The spectrum of physical fluctuations of these fields is dual to the 
spectrum of glueballs in the gauge theory, as these fields come from the closed sting 
sector. The physical massive fluctuations of the minimal metric+dilaton system 
consists of one spin-2 mode (5 degrees of freedom), and one spin-0 mode.^'' The 
fluctuations of the axion field correspond to pseudoscalar glueballs. They do not 
mix with those in the scalar sector of the metric-dilaton system, since we neglect the 
backreaction of the axion on the geometry^^. 

Throughout this section we consider only the IR asymptotics of the type (|3.13| ), 

^(r)~-(^)%..., a>l (6.24) 

with no further assumptions on the subleading behavior. We will consider the case 
with singularity at finite r in Section W. 



6.5.1 Scalar glueballs 

In 5D Einstein-Dilaton gravity there exists a single gauge invariant spin-0 mode^^, 
C(r, x), satisfying the equation (|6.2|) with 

Eo(r) = ^-A{r) + i \ogX\ M{r) = 0, (6.26) 



"'^'^See e.g. ||Tl| for a complete discussion of the identification of the physical fluctuations and the 
corresponding field equations. In the massless sector there are a massless spin-2 (2 polarization), 
one massless spin-1 (2 polarizations) and 2 massless spin-0 modes. However we will not have 
massless modes in our spectra, so we will not consider this case further. In [ pl| it was shown that 
in general the presence of a massless spin-2 mode is only possible if an IR singularity appears and 
if special non-local boundary conditions are put at the singularity. This is compatible with the 
Weinberg- Witten theorem Q . 

"'^^However, they are expected to mix with 77' if we introduce flavor branes. 

-'^^Here "gauge invariace" refers to the linearized 5D diffeomorphisms. The precise definition of 
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The effective Schrodinger potential is given by eq. (|6.6|) . Notice that, both for 
large and small r, the second term in Bq ( |6.26|) is negligible. Therefore the leading 
asymptotics are, 



We have a mass gap and discrete spectrum if and only if a > 1. 

In the UV, the gauge invariance of ( indicates that it is dual to the renormaliza- 
tion group invariant operator /3(A)Tr[F^] @]. 



6.5.2 Tensor glueballs 

The massive spin-2 glueballs are described by transverse traceless tensor fluctuations 
hij of the 4D part of the metric: 

ds^ = e2^(") {dr^ + {7]ij + hij)dx'dx^) (6.28) 

These fluctuations satisfy the equation (|6.2|) with 



B2{r) = ^A{r), M{r) = (6.29) 
The effective Schrodinger potential has the same asymptotics as (|6.27|) : 



9 
4 

V 

4 f.2 ) 



V2{r)-{f;\ ^ (6.30) 



Together with (|6.27|) this confirms (|6.22|) . However, due to the difference between 



( |6.26|) and ( |6.29| ), the spin-0 and spin-2 glueball spectra are not degenerate. This is 
unlike the AdS/QCD models with exact AdS metric and constant dilaton: there the 
scalar and tensor modes are exactly degenerate. We will see in an explicit background 
that the lowest-lying spin-0 glueball is lighter than the lowest spin-2 glueball. We 
expect this fact to be generic, although we can not provide a proof in our set-up. 

6.5.3 Pseudo-scalar glueballs 

The Einstein frame axion action in the conformal coordinates reads: 

Sa = — — / d''xZ{X)e^^{daf. (6.31) 



this field is 

^ = ^ - WtT^'^ = ^ - fl^^"^' (6-25) 
3X{r) /3(A) 

where Sfjj and ip are the fluctuations in the dilaton and in the scalar part of the gij metric component. 
See e.g. p. 
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Since the axion appears quadratically, this is also the action for the fluctuations. We 
thus have: 

B,{r) = -A{r) + -\ogZ{X). (6.32) 

To leading order in string perturbation theory, Z{\) = X^. However, this in general 
is expected to receive corrections from the 5-form, similar to the dilaton potential. 
Indeed, if this were not the case one would find a puzzling result: one would obtain 
a continuous spectrum for the pseudo-scalar glueballs starting at m = 0. To see this, 
assume as in Section | that Z{\) = X'^ for large A. Then, using eq. ( p.l5| ) in ( |6.32| ) 
we obtain: 

Ba{r) =l(l-^)A{r) + ^^-{a-l)\ogr/R 



{d - 2) {r/Rf d^2 



Ucy-l)\ogr/R d = 2, ^^-^^^ 



where we used ( |6.24D . The IR asymptotics of the Schrodinger potential are (using 

V(r)^lrAid-^r/R']ir/Rf-'^ d^2 

U^(«-l)^-!(«-l)]^ d = 2. ^'-'"^ 
Thus the potential and the spectrum have the same features as the other glueballs, 
unless the perturbative result d = 2 is unmodified. 

The asymptotic mass ratio for large n of the to 0++ glueball states can be 
read-off comparing the large r asymptotics of (|6.34|) for d ^ 2 and (|6.27|) : 

Using the expected asymptotic glueball universality argument (as in ( |6.22|) ) 

<^oo(0++) <_(2++) ' ^ • ^ 

we can determine 

d = 4. (6.37) 
This result predicts an interesting renormalization of the bare axion kinetic term. 



( lOT 



It is appropriate to point out that the effective Schrodinger potential for the 
trajectory of glueballs can be written in terms of the background axion solution 
( CT) as 

K(r) = \- (6.38) 
4 a 

An interesting corollary of this relation is that the potential is independent of the 
UV e-angle of QCD, Ouv 
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When 9 ^ Q glueballs of different parities mix as well as their energies become 
^-dependent to next order in 1/iVc. The 6'-dependence of the glueball spectrum can 
be calculated by considering the first order ((9(l/iV^)) backreaction of the axion 
solution to the QCD vacuum. This can be an interesting test as lattice data on this 
exist, |2^. We will not attempt however the calculation in this paper. 



6.6 Singularity at finite ro 

In the previous subsections we considered backgrounds with infinite range in r. Here 
we discuss the case in which the IR singularity is at some finite r = tq. As discussed 
in Appendix ^ and summarized in Table these backgrounds generically lead to a 
confining string potential. To analyze the mass spectrum, consider the case when 
the IR singularity has the following form: 

A(r) ~ (5 log(ro — r), r ^ Tq. (6.39) 

The effective Schrodinger potential ( |6.6| ) has the same asymptotic form in the IR 
both for the scalar and the tensor glueballs. This is because the functions (6.26) and 



( |6.29| ), differ only by a function of X{r) which, as shown in Appendix asymptotes 
to a ((5-dependent) constant as r — > tq. Then, both for the spin-0 and the spin-2 
glueballs, the effective Schrodinger potential has the following asymptotic form: 

l-(r)~^l (r-0), V{r)^\^^^^, (r ^ oo). (6.40) 
4 4 (r — ro) 

For 5 > 2/3, V +00 in the IR, and by the same general argument we used in 
subsection 4.1 we obtain a mass gap and a discrete spectrum. The treatment of the 
case < 5 < 2/3 (in fact < 5 < 1) requires extra care, as we discuss in the next 
subsection. 

The large mass asymptotics of both the scalar and the tensor glueballs in the 
backgrounds (|6.39|) are universal. They depend neither on 5 nor the details of the 



metric in the bulk: due to ( |6.40| ), the Schrodinger equation for large eigenvalues is 
effectively the one for a particle in a box of size tq, so for large mass eigenstates we 
obtain 



n 



2 



< ~ -• (6.41) 

This does not a priori prevent the mesons to have a linear mass spectrum, since this 
is driven by the tachyon dynamics^^, as in the infinite range case. 
In the case of power-law behavior for A(r), 

C 

A(r) 6,C>0 (6.42) 



^"This observation avoids the arguments put forward in regarding the meson spectra in 
gravity duals. 
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the potential in the IR always asymptotes to +00, and it is steeper than [vq — r) ^ 
in the IR: 

6.6.1 The pathologies for < 6 < 1 

As discussed in [|ll|] in a different context (see also |^6| for a related discussion), this 



range of parameters is somewhat pathological, since it requires additional boundary 
conditions at the singularity, and the spectrum is not determined by the normaliza- 
tion condition alone. 

The Schrodinger equation for a generic mass eigenstate close to tq is: 



-^ + V(r) ^ -^ + ^^ — ^^ = m2^. (6.44) 
4 (r — ro) 



For r ~ To we can neglect the mass term on the r.h.s, and find the asymptotic solution 
close to Tq: 

Hr) ~ Ci(ro - r)3^/2 + C2(ro - r^-^^^^ (6.45) 

For 6 < 1 both solutions are square-integrable, and they both vanish at Tq if in 
addition 6 < 2/3. Therefore, for < 5 < 1, normalizability alone is not enough to 
fix the spectrum uniquely. One has to specify some extra boundary conditions at 
the singularity, which may be given by fixing the ratio Ci/c2.^^ In contrast, for 6 > 1 
normalizability in the IR forces the choice C2 = 0, and there is no ambiguity. 

Ultimately it is this extra input at the singularity that determines the spectrum 
in a background with 6 < 1, and not the dynamics of the theory at any finite energy. 
This situation is not so different from the hard- wall models |0, Q , where one also has 
to specify IR boundary conditions for the fluctuations to compute the spectrum. 

We note here that the background studied of Csaki and Reece in [O] falls in this 



class of examples: one can easily check that its metric in conformal frame behaves 
as in eq. ( |6.39|) with 5 = 1/3. In computing the spectrum, IR Neumann boundary 
conditions are chosen in [jl5l , but according to the present discussion this is as good 
a choice as any other. 

6.7 Adding flavor 

A small number Nf -C Nc of quark flavors can be included in our setup by adding 
space-time flUing "flavor-branes" . In this case they are pairs of space-fllling DA — DA 
branes. It was proposed in that the proper treatment of the flavor sector (including 
chiral symmetry breaking) involves the dynamics of the open string tachyon of the 
DA — DA system. According to this, the meson sector of the 4D gauge theory is 



^^In operator language, the Hamiltonian of this problem is symmetric but not essentially self- 
adjoint, and it admits an infinite number of self-adjoint extensions, each with a different spectrum, 
parametrized by the choice of Ci/c2- 
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captured holographically by the open string DBI+WZ action, which schematically 
reads, in the string frame. 



S[T,A'',A''] = Sdbi + Swz (6.46) 
where the DBI action for the pair is 

Sdbi = j drd'x y Str V{T) det {g^, + D^^TW^^T + F^,) + (6.47) 

+ det [g,, + D{,TW,yT + F« )^ 

Here T is the tachyon, a complex Nf x Nf matrix. Aj^'^ are the world-volume 
gauge fields of the U{Nf)L x U{Nf)R flavor symmetry, under which the tachyon 
is transforming as the {Nf,Nf), a fact reflected in the presence of the covariant 
derivatives^^ 

D^T = d^T -zTAf^ + iA^T , D^^T^ = d^T^ - lA^^T^ + iT^ A^ (6.48) 

transforming covariantly under 

T^VnTVl , A"^ ^VL{A'^-tVldVL)Vl , A« ^ r^(A^ - 2l^^dF«)V^^ (6.49) 

as well as the field strengths F^'^ = dA^ ji — iAl^r A A^ ji of the A^'^ gauge fields. 
A = e* = Nce'^ is as usual the 't Hooft coupling. We have also used the symmetric 
trace (= Str) prescription although higher order terms of the non-abelian DBI action 
are not known. It turns out that such a prescription is not relevant for the vacuum 
structure in the meson sector (as determined by the classical solution of the tachyon) 
neither for the mass spectrum. The reason is that we may treat the light quark masses 
as equal to the first approximation and then in the vacuum, T = rl with r real, and 
this is insensitive to non-abelian ramifications. Expanding around this solution, the 
non-abelian ambiguities in the higher order terms do not enter at quadratic order. 
Therefore, for the spectrum we might as well replace Str —>■ Tr. 
The WZ action on the other hand is given by^^: 



Swz = T4 C A Str exp [^27^Q;'JF] 



(6.50) 



where M5 is the world- volume of the D4 -D4 branes that coincides with the full space- 
time. Here, C is a formal sum of the RR potentials C = X]n(~^)'^^™' 
the curvature of a superconnection A. In terms of the tachyon field matrix T and 



We are using the conventions of §. 

This expression was proposed in and proved in ^] using boundary string field theory 
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the gauge fields and living respectively on the branes and antibranes, they 
are (We will set 27ra' = 1 and use the notation of |1T8[| ): 

The curvature of the superconnection is defined as: 

J^ = dA-iAAA , dj^ -iAAJ^ + iJ^ AA = (6.52) 
Note that under (flavor) gauge transformation it transforms homogeneously 

In the relevant definitions and properties of this supermatrix formalism can be 
found. 

By expanding we obtain 

Swz = T4 y" C5 A Zo + Cs A ^2 + Ci A Z4 + C_i A Ze (6.54) 

where are appropriate forms coming from the expansion of the exponential of the 
superconnection. In particular, Zq = 0, signaling the global cancellation of 4-brane 
charge, which is equivalent to the cancelation of the gauge anomaly in QCD. Further, 
as was shown in 

Z2 = d^i , ni = tSTr{V{T^T))Tr{AL - A^) - \ogdet{T)d{StrV{T^T)) (6.55) 

This terms provides the Stuckelberg mixing between Tr[Aj - A^] and the QCD 
axion that is dual to C3. Dualizing the full action we obtain 

Scp-odd = ^ J d'x,/^Z{X) {da + i^lif (6.56) 
= ^ / d'x^gZ{\) (^,a + Cd,V{r) - ^I^ViT)A^'j 

with 

C = 5>logdetT , AL-An^^AH + {Al-A%)X'^ (6.57) 

and where we have set the tachyon to it vev T = rl . This term is invariant under 
the U{1)a transformations 

C-C + e , A^^A^~J^^d,e , a^a-NjeVir) (6.58) 
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reflecting the QCD U{1)a anomaly. It is this Stuckelberg term together with the 
kinetic term of the tachyon field that is responsible for the mixing between the 
QCD axion and the 1]'. In terms of degrees of freedom, we have two scalars a, ( 
and an (axial) vector, A^. We can use gauge invariance to remove the longitudinal 
components of A^. Then an appropriate linear combination of the two scalars will 
become the O""*" glueball field while the other will be the rj'. The transverse (5d) 
vector will provide the tower of U{1)a vector mesons. 

The next term in the WZ expansion couples the baryon density to a one-form 
RR field Ci. There is no known operator expected to be dual to this bulk form. 
However its presence and coupling to baryon density can be understood as follows. 
Before decoupling the Nc -D3 branes, its dual form C2 couples to the U(l) b on the 
D3 branes via the standard C2 A Fb WZ coupling. This is dual to a free field, the 
doubleton, living only at the boundary of the bulk. Once we add the probe + D4 
branes the free field is now a linear combination of A^ and an Nf/N^ admixture 
of A^ originating on the flavor branes. The orthogonal combination is the baryon 
number current on the flavor branes and it naturally couples to Ci. Therefore the 
Ci field is expected to be dual to the topological baryon current at the boundary. 

Finally the form of the last term requires some explanation. By writing Zq = dfl^ 
we may rewrite this term as 

y Fo A ^5 , Fo = dC_i (6.59) 

Fo ~ Nc is nothing else but the dual of the five-form field strength. This term 
then provides the correct Chern-Simons form that reproduces the flavor anomalies of 
QCD. Its explicit form in terms of the gauge fields Al b and the tachyon was given 
in equation (3.13) in 

To proceed further and analyze the vacuum solution we set T = r 1 and set the 
vectors to zero. Then the action (|6.46|) collapses to 

S[r,AM] = NcNf J drd\ e-^'V^r)^ - det {g^, + d^rd.r) (6.60) 

Following we assume the following tachyon potential, motivated/calculated in 
studies of tachyon condensation: 

V{t) = Voe-'^^' (6.61) 

where fi has dimension of mass. It is fixed by the requirement that r has the correct 
bulk mass to couple to the quark bilinear operator on the boundary. 

In our minimal setup, the brane-antibrane system fills the whole bulk. Therefore 
these fields are bulk fields. We will eventually expand the action at most to quadratic 
order in the gauge fields. 
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Chiral symmetry breaking in the IR is described by a non-trivial tachyon profile. 
For small Nf we can neglect the backreaction of the tachyon on the metric-dilaton 
system, and solve the equation for the tachyon profile on a given background, e.g. 
one of the confining backgrounds we discussed. Once a solution for the tachyon 
is found, the spectrum of mesons is given by the spectrum of fluctuations around 
this background. For example, vector mesons are described by the fluctuations of 
the components Ai around the Ai = configuration, in a given background for the 
metric, dilaton and tachyon. 

6.7.1 Tachyon dynamics 

In the conformal frame, the action (|6.6CI| ) becomes: 



S[t] = KNf j rfrdW^^('^)-*('-V(r)-^e2^^M + f{r)\ (6.62) 

from which we obtain the nonlinear field equation: 

f + (?,As - f + e^^^V + AAs - $ {ff + /iV(f)2 = 0. (6.63) 

The tachyon is dual to the dimension 3 quark bilinear operator. Near the boundary, 
r — > 0, we expect r = mr + crr^ + . . Thus, in the UV We may ignore the non- 
linear terms in eq. (|6.63| ), which then reduces to the equation for a free scalar field 
with mass /i on an asymptotically AdS^ background. In order for this to be dual 
to the quark bilinear operator, with naive dimension 3 (to leading order), we need 
3 = 2+ - i^H'^, hence fiH^ = 3. 

It is argued in P that consistency of the bulk gauge theory (i.e. absence of 
extra gauge anomalies in the IR ) requires the tachyon to diverge before or at the 
singularity. In Appendix we analyze the possible singularities of the solutions 
of eq. (|6.63| ), in backgrounds with IR asymptotics (|6.24|) . We show that the only 
consistent solution for r — >■ oo, is such that the tachyon diverges exponentially: 



r(r) ~ To exp 



2R 



oo, (6.64) 



where tq is an integration constant determined by UV initial conditions. 

We also analyze possible singularities of the solutions at finite r. We find that 
generically, the tachyon cannot diverge at any finite r, where both As and $ are 
regular, except special points where iAg — $ = 0. This does not happen in our 
backgrounds. Instead, the generic solution of ( 3.63| ) has a singularity at finite r*. 



where T(r*) stays finite but its derivatives diverges: 

r ~ + 7^^"* - r. (6.65) 



^■^For simplicity, here we take all quark masses to be equal. The tachyon is therefore proportional 
to the identity matrix in flavor space. 
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Such solutions are unphysical, since around r* the backreaction on the metric is no 
longer negligible: the tachyon stress tensor diverges as l/(r,, — r), and our assumption 
that the tachyon does not perturb the background is invalid. On the other hand this 
is not physically reasonable, since adding a small number of flavors should not change 
dramatically the pure gauge dynamics in the large Nc limit. 

Discarding all but the exponentially divergent solution singles out special initial 
conditions in the UV, which correspond to fixing the chiral condensate as a function 
of the quark mass f^, i.e. the coefficients of the subleading and leading terms in the 
UV expansion of T(r). 

6.7.2 Vector mesons 



Once the correct tachyon profile is found from eq. ( |6.63| ), this enters the action for 



the tachyon and the bulk gauge fields fluctuations, and determines their spectrum. 
The resulting 4D mass eigenstates correspond to the various mesons in the dual 
theory. Here, we only consider the vector mesons, that correspond to the transverse 
vector components of the 5D gauge fields, Ai = + Af. 

The quadratic action for the gauge fields is, from eq. (|6.46|) : 



S--\j drd'x e-''V{T)y^r''9'''F^,F,^, (6.66) 

where g is the effective (open string) metric felt by the gauge fields in the presence 
of the tachyon: 

ds^ = (e^^s + (f)2) dr^ + e^^''r]ijdx'dx^. (6.67) 

This metric is still asymptotically AdS, since e'^^'^ dominates in the UV, however, 
although still conformally fiat, it is not in the conformal frame. It differs considerably 
from the bulk background metric in the IR. 

The large r behavior of As{r) and r(r) are, from eq. ( ^.161 ) and ( |6.64| ): 



a - 1 

As{r) — - — \ogr/R, r(r)~roexp 



2R 



a>l. (6.68) 



The second term dominates grr in the infrared. To recast the action in the form (|6.1| ), 
and read-off the effective Schrodinger potential for the mesons, we change variables 
from r to r. Using ( |6.68| ) to express As as a function of r in the IR, the effective 



metric becomes for large r: 



ds"^ ~ (ir^ + ( ^^log^/^o ) r]ijdx'dx^, (6.69) 



^^Notice that the backreaction is not problematic if the tachyon itself, and not just its derivative, 
diverge: the stress tensor is multiplied by the tachyon potential, that vanishes exponentially fast 
as T — > c», resulting in the recombination of the branes-antibrane pairs in the IR, which leaves the 
unperturbed metric and dilaton background. 
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where we have neglected the first term in grr- We now pass to a new conformal 
frame, by changing variables from r to f, defined by 



a£2 



(a-l)/2 



dr= [^logr/To] dr + ■ ■ ■ , (6.70) 



which is solved asymptotically for large r by: 

+ ■■■. (6.71) 



2_^2\ (a-l)/2 



(logr/ro)("-^^/^ 

To leading order we can also replace logr/ro by log t/tq in the above relation and 
the metric reads: 



ds^ = e^^^''^ [dr^ + riijdx'dx^) 



■logr/ro 



The action for the transverse vector fluctuations becomes: 



(rff2 + r]ijdx'dx^) . (6.72) 



S = -\j dfd^xe-''V{f)e^^'^ [{dfAi)^ + (9, A)'] , (6.73) 



and has the same form as in (B.ll) with 



^ -^(^) ^ '^(^) + 1 log V{r{r)) (6.74) 

Asymptotically the last term dominates (it behaves like r^, which is exponential in 
the original r coordinate, while As grows logarithmically and $ a power-law of r), 
and we find, using eq. (|6.61|) : 

B{f) ~ - A (^^y ' (logf/ro)--^ (6.75) 

/^From the general analysis of section |6.1| , and in particular from eq. ( p.6|) , the lead- 
ing behavior of the vector meson Schrodinger potential is that of a (logarithmically 
corrected) harmonic oscillator, therefore it exhibits an approximately linear mass 
spectrum^^. This is a concrete realization of the general mechanism described in [^]. 

Notice that the meson spectrum is generically controlled by a different energy 
scale than the one that sets the glueball masses: the two scales are 

-1 _ 3 / «£^y"-'^/' i/^iV"' 

■^'■glueballs j-^j -^'■mesons ^ I ) ^ R \ R ) ' 1^0. lOj 

Interestingly, the two scales happen to coincide in the special case a = 2, in 
which the asymptotic glueball spectrum is also linear. 

As a final remark we comment on the importance of corrections due to 

the bulk (closed sector) on the meson spectrum. It seems that at least for questions 
of the meson spectrum, lattice calculations indicate that such corrections are small. 
The errors for the pion and p-meson masses were estimated at around 4% for A'^c = 3 
in 



|6[ working in the quenched approximation. 



26 



One can get rid of the extra log by a slight modification of the tachyon potential. 
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7. The parameters of the correspondence 



QCD with gauge group SU{Nc) has three parameters: the bare couphng constant 
Aq, the theta-angle 6 and the number of colors N^. On the other hand, through 
dimensional transmutation the bare coupling constant is replaced by the dynamically 
generated strong coupling scale Aqcd- In the large Nc limit therefore, apart from 6, 
this is the single parameter of the theory. It sets the scale for the glueball and meson 
masses. 

7.1 Parameters in the gravitational action 

On the gravity side, we have a number of parameters entering the Lagrangian (and 
the fundamental string action) and a set of integration constants for the vacuum 
solutions. The parameters that appear in the action are the 5D Planck scale M, the 
string scale 4, the AdS radius i (via the overall scale of the potential, Vq = 12/£^). 
Moreover, the potential as a function has dimensionless parameters. In its weak 
coupling expansion, they are in one to one correspondence with the coefficients of the 
/3-function, 6„. Among the parameters 6„, as we discussed above, the only scheme 
independent coefficients are bo and bi. For practical reasons we restrict ourselves 
to potentials parameterized only by these two. Moreover, the ambiguity in the 
normalization of A, amounts to bo being a parameter that we eventually fit. bi/b^ 
however remains and we take its value from the QCD /3-function. The strategy 
however is clear: although the potential is an a priori arbitrary function, our eventual 
choice will have very few adjustable parameters. 

The Planck scale governs the strength of bulk interactionss. In the underlying 
string theory, it is determined in terms of the string scale. However here, this relation 
is not known and it will have to be taken as an extra parameter, that can be fixed 
by matching to interactions, or to the finite temperature free energy. 

On the other hand, either of the two dimensionfull parameters i or is can be 
used to set the units. We choose to measure all the dimensionfull quantities in units 



Apart from the parameters of the action, there are in general three integration 
constants that parameterize the solutions to the Einstein-dilaton system. 

Now, we discuss them one by one. In the dual gauge theory the /3-function is 
fixed, and this is in one-to-one correspondence with the function X or the superpo- 
tential W on the gravity side. Once the potential V is given, there are infinite number 
of W that solve ( |2.10|) , parameterized by a single boundary condition. However in 
appendix (^ we show that the confining asymptotics of the form, 



of i. 



W 
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Figure 3: The profile of tlie scale factor and the dilaton in a typical solution. is the 
curvature singularity where the scale factor shrinks to zero and the dilaton blows up. 

is a unique solution to the equation ( |2.10| ). Therefore the requirement of confine- 
ment uniquely fixes the superpotential W and reduces the number of independent 
integration constants from three to two. 

The IR asymptotics above indeed evolves into an asymptotic AdS space in the 
UV as, 

Vo^ + 0(A). (7.2) 
Given the superpotential, the equations of motion reduce to two first-order equations 

m- 

7.2 Reparametrization symmetry and integration constants 

The remaining two initial conditions of the first order system of motion are related 
to the "radial reparametrization symmetry" of the equations of motion. Indeed from 
the Einstein's equations ( |2.9[ ) one learns that, for any solution A^{u), $=„(u) there 
exist other solutions parametrized by two numbers: 

A{u) = A^{u - uo) + Ao, $(«) = $*(m -Mo)- (7.3) 

The parameters Aq and Uq are in one-to-one correspondence with the two integration 
constants. 

It is useful to describe the symmetries ( [7.3| ) graphically. In fig. |^ we exhibit a 
typical solution that obeys our UV and IR criteria. Given the solution in fig. |^, one 
can generate another by shifting A{u) vertically. This corresponds to the first shift 
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symmetry A A + Aq in {7_^) . The second shift symmetry u ^ u + uq, corresponds 
to generating another solution by shifting both X{u) and A{u) horizontally in fig. ^ 
In practice, the solution is fixed by two initial conditions set at an arbitrary point 
Uq: {A{uo),X{uo)} = {Ao,Xo}. 

The next observation is that the integration constant uq is just a gauge artifact. 
This is because the original system is translation invariant in the radial variable. 
In fact, using this reparametrization symmetry, one can use A in place of the radial 
variable u and then the physics is completely specified by the function A{X). Clearly a 
horizontal shift in fig. ^ although changes Uq, leaves A{X) invariant. In other words, 
the only physical integration constant in the system is the one that parametrizes the 
solutions of 

dA/dX = (3~\X). (7.4) 



The general solution is: 



and it is invariant under: 



Xo^Xo + 5, Ao^Ao+ I -j^. (7.6) 



Ao 

Therefore, in (|7.5| ) there is a single combination of Aq and Aq that specifies the full 
solution. This can be taken as the value of A at some fixed A, e.g. Aq = A{Xo). 

Therefore, all physically distinct solutions only differ by a constant shift in A. 
Fixing the integration constant Aq = A{Xq) is equivalent to specifying Aqcd in the 
gauge theory, because it sets the energy scale through the relation E = expA. In 
fact, this is the only way the integration constant Aq affects any physical quantity: 
a change in Aq induces a constant rescaling of all dimensionfuU quantities, such 
as masses, confining string tension, etc. In particular, mass ratios are completely 
independent of all integration constants, and only depend on the parameters that 
appear in the gravity Lagrangian. 

Consider for example the scalar or tensor fluctuations. The corresponding spec- 
tral equations follow from an action of the form (|6.1|) , with the functions B and M 
given by eqs. ( |6.26|) and (|6.29| ). We can change coordinates from r to $ = log A in 
(|6.1|), and derive the corresponding spectral equation for the fluctuations C(^)- 

—2A 

-e-^^'^d^ [e'^^^d^Wd^O = (7.7) 

Under a constant shift A A + 5Aq, and B also shifts by a constant. The left hand 
side is therefore invariant, and the right hand side is rescaled by e~'^^'^°\ Thus, the 
only effect on the spectrum is an overall rescaling of all the mass eigenvalues by e^^° . 



^^Note that this is not the case for a conformally invariant theory. In that case, one cannot use 
the couphng as a new coordinate as it is constant. 
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The same considerations hold for the confining string tension, eq. (|3.12 ): a 



constant shift in A does not change the position of the minimum of As{r), but it 
only rescales the tension by 6"^^^°. In particular, the ratios m^/Tg are independent 
of the integration constant. 

A change in Aq has the same effect also on the scale governing the perturba- 
tive running of the coupling, that gives rise to dimensional transmutation: if initial 
conditions are chosen at Aq ^ 1, integration of the /9-function equation ( [7.5|) leads 
to: 

11 E e^o 

We can identify the "perturbative" RG-invariant QCD scale as follows. Integrat- 
ing eq. (|7.5| ) for small A, Aq up to two loops, with /5(A) ~ — feo-^^ ~ ^i-^^ + Oi^X^), we 
obtain: 

^(A)-7^-^log(&oA)+0(A) = A^--^ — ^log(6oAo) + 0(Ao) - constant. (7.9) 



Therefore, the scale 



Ap = -- 



lexp 



(7.10) 



(feoAo)"^/"" 

is approximately independent of Aq as long as it is small. This scale appears in the 
UV expansion of the coupling in the form: 



1 E E 

-TT = logT 72 log log T 

boX Ap ¥q Ap 



log 7 l2 loglog— + • • • (''•ll) 



It is the same scale appearing in the UV expansion of the solution in conformal 
coordinates, (|2.14|) , as one can see by substituting E ~ 1/r on the l.h.s. of eq. 



All the different scales we have analyzed above behave in the same way under a 
change in the integration constants, so the relations between them is a property of 
the gravity model, not of each particular solution. 

In the explicit examples we present in Section |^, we fix the energy scale to match 
the lowest glueball mass, and as we have discussed this fixes unambiguously all other 
dimensionfull quantities. In particular we obtain a value for Ap in eq. ( [7.10|) : 

Ap = 290MeV. (7.12) 

This value is larger than the usual QCD value (~ 200 MeV). However one should 
keep in mind that the definition of the strong coupling scale in perturbation theory 
is somewhat arbitrary, and moreover we are not including the effect of quarks in 
the running of the coupling. When comparing with data (lattice or experiment) it 
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is more meaningful to look at unambiguous quantities, e.g. the value of the strong 
coupling constant as at a given energy. For example, we find^^. 

a,(1.2Ge1/) = 0.34 (7.13) 

which is very close to the experimental value a's^^\l.2GeV) = 0.35 ± 0.02 

In summary, like in QCD, in the gravity side solution is specified by the (3- 
function plus a single dimensionfull quantity e^°/i, that parametrizes the different 
solutions and sets all the relevant mass scales. It can be related to Aqcd of the gauge 
theory, and it can be understood holographically as the constant of motion that is 
preserved under the shift symmetry u ^ u + uq. This is precisely Aq in our set-up. 



8. Concrete backgrounds 

In this section we present explicit backgrounds that exhibit all of the features we 
require (asymptotic freedom, confinement, discrete spectrum). Then we compute 
the glueball spectra numerically. 

We consider two backgrounds belonging to two distinct classes. The first is a 
background with an exponentially decaying scale factor, and with an infinite range 



of the conformal coordinate. We focus on the case a = 2 in (see equation ( 6.24 )). As 
shown in section |6.3| this gives an asymptotically linear glueball spectrum. Secondly, 
we analyze an example of a background with finite range of the conformal coordinate. 
In both cases we fix the 5D theory by providing a function X{X) that interpolates 
between the required UV and IR asymptotics. As we discussed, this is equivalent 
to fixing the exact /9-function. The RG-fiow trajectory is further specified by the 
UV initial conditions, which we input for the numerical integration. This fixes the 
gravity dual completely. 

In this paper we only present the glueball spectra. Although straightforward 
in principle, the meson spectra require considerably more complicated numerics. 
The main obstacle from the numerical point of view is identifying the correct initial 
conditions for the nonlinear tachyon equation, (which is then used as an input in the 
computation of the meson spectrum). Therefore we leave the computation of the 
meson spectrum for future work. 

Finally, we compare the glueball spectra with the available lattice data. For the 
model with infinite range of r and a = 2 we can fix the parameters in such a way to 
produce a very good agreement, at a quantitative level. 



^^See section ^ for a more detailed discussion, and in particular for the relation between A and 

a. 
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8.1 Background I: unbounded conformal coordinate 

For an asymptotically free, confining theory, the function X{X) has the following UV 
asymptotics (see eq. ( |2.18| )) 

X(A) ~ -yA- ^A2 + ... A^O (8.1) 

where bk are the k-th order coefficients of the perturbative /3-function. In the IR we 
require (see ( |3.17|) ): 

where the parameter a determines the large-r behavior of the scale factor: 

A~-Cr" a = -^^. (8.3) 

8 a 



We seek for a function of A that interpolates between the two asymptotics ( p.l|) and 
( p.2|) . A simple function that is regular and has this property is, 

X(X) = —^ (2^0 + 3^0^' (SA) 

^ ^ 3 + 2boX 9(1 + A2) (1 + gL (262 + sb,) log(l + A^)) ' ^ ' ^ 



This expression is motivated by the UV and the IR asymptotics in (|8.1| ) and 
^ and by the requirement that there are no poles or branch cut singularities in 
A. Also, the function X(A) (hence also /5(A)) is strictly negative for A > 0, therefore 
there are no IR fixed points. 

Starting from eq. ( p.4|) , we solve for the metric and dilaton using eqs. ( |2.12|) : 

X = -±X{X)W{X)Xe^, i = -liy(A)e^. (8.5) 
The superpotential W{X), is given in terms of X as in (|2.22| ), 



2a/ 3 



(8.6) 



and in writing ( ^.5] ) we have explicitly extracted the overall scale £. In the integration 
of (|8.5| ), we fix the integration constants as: 

A{ri„) = Ao, X{rin) = Aq. (8.7) 

for Tin/i ^ 1 and Aq ^ 1, in order to implement the correct UV asymptotics. 

The scalar and tensor glueball spectra are completely fixed by the metric and 
dilaton background. For the pseudoscalar glueballs, we need to specify also the axion 
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Figure 4: The scale factor and 't Hooft coupling tiiat follow from ( |8.4D , bo 
initial conditions Aq = 0, Aq = 0.05 at r = 0.36. The units are such that i 
dashed line represents the scale factor for pure AdS. 



= 4.2, and 
= 1. The 



kinetic function Z{^) appearing in eq. (|6.31 ). From the discussion Sections 5 and 
6.2.3, we assume the asymptotic behaviors: 



We take the function Z{\) to be the simplest one satisfying these asymptotics 

Z{\) = l+CaX^ (8.9) 

where we have fixed the = 1 by an overall rescaling of Z{X), which does not affect 
the glueball spectrum, is an extra dimensionless parameter^^. 

8.1.1 The glueball spectra in background I 

We solve the eq. ( |6.5D with the Schrodinger potential ( |6.6|) numerically. We compute 
the spectrum of scalar and tensor glueballs where the function B in ( |6.6|) is given by 
eqs. (|6.26| ) and (|6.29|) respectively, whereas the 5D mass-term M in ( |6.6| ) is zero. 



One has to supply the Schrodinger equation with the boundary condition in the 
UV, (as r ^ 0), 

^^Cori + Cir-i (8.10) 

Particle states correspond to normalizable solutions. Therefore, normalizability in 
the UV requires Ci = 0. Normalizability in the IR, on the other hand fixes the 



^^We could take a more general form that includes the pcrturbative string theory term ~ A^, 

Z{X) = l+baX^ +CaX\ 

however for the sake of simphcity we set 6a = in our fits. A non-zero ba would imply a different 
preferred value for Ca, but this does not change the spectrum significatively. However this could 
have a non-negligible effect on the axion profile, see Fig. |l|. 



discrete values for m in (3^). In practice, we use the shooting method to determine 
the spectrum: we scan the values for m and pick the values at which an extra node in 
the wave function appears. Precisely at this value of m, the wave function becomes 
normalizable in the IR. 

In principle, the spectrum depends on the parameters of the background, Bq and 
61, the integration constants of the geometry Aq and Aq (eqs. (|8^ )) and the boundary 
condition of ( ^.10|) , i.e. Cq. However, not all of these parameters affect the spectrum 
nontrivially. 

• The constant Cq is clearly immaterial, due to the linearity of the equation for 
iplr). We set Co = 1 without loss of generality, 

• As we discussed in Section 7, the only physical integration constant for the 
background is the choice of Aq at some value Aq, and it only affects the overall 
scale of the masses. This expectation is confirmed by the numerical results, as 
shown explicitly in figure ^. Thus, the mass ratios will be independent of Aq 
and Aq, as well as of the AdS scale. The overall energy scale can then be fixed 
by matching e.g. the mass of the lowest state in the spectrum. 



"'0 
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(a) (b) 

Figure 5: Dependences on initial condition Aq of (a) the mass ratios Rqq = ?n,o*++/mo++ 
(squares) and R20 = 77i2++/mo++ (triangles); (b) the absolute scale of the lowest lying 
scalar glueball (shown in Logarithmic scale). 



As discussed in [||], 60 cannot be determined from first principles in our setup, 
as the overall coefficient in the relation (|2.16 ) between the dilaton and 't Hooft 
coupling is not known. We keep bo as a free parameter. On the other hand, 
the ratio is independent of such normalization. In pure YM this ratio is 
given by 51/121 and this is what we use. 
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The superpotential completely fixes the scalar and tensor glueball mass ratios; 
the pseudoscalar glueball masses also depend on the additional parameter Ca 
that enter the definition of the axion kinetic term, eq. ( |8.9| ) 



In light of the above, we vary only fog for the purpose of fitting the scalar and 
tensor glueball lattice data, and Ca to fit the pseudoscalar glueball data. 



V[r] 



. 




Figure 6: Effective Schrodinger potentials for scalar (solid line) and tensor (dashed line) 
glueballs. The units are chosen such that i = 1. 




ID 2D 30 40 5D 60 70 2 4 6 



(a) (b) 

Figure 7: (a) Linear pattern in the spectrum for the first 40 O"*"^ glueball states. is 
shown units of 0.007^"^. (b) The first 8 0++ (squares) and the 2++ (t riangles) glueballs. 
These spectra are obtained in the background I with bo = 4.2. 

We perform most of the numerical analysis for the background that gives linear 
spectrum, i.e. a = 2 (we discuss the dependence of the spectrum on the parameter 
a at the end of this section.). To make the numerics easier, we fix £ = 1 and work 
in dimensionless units. The geometry looks typically like in Figure ^, the effective 
Schrodinger potentials as in Figure ||, and the glueball spectrum as in Figure |^. 

We note that, unlike the simple AdS/QCD setup, the scalar and tensor glueballs 
are not degenerate, but the tensor glueballs are generically heavier than the scalar 
ones with the same quantum number n. The tensor-scalar mass difference decreases 
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for larger n, indicating that the slopes governing the asymptotics of the two spectra 
are the same. This is in accord with our discussion in section 



8.1.2 Lattice Data 



Available sources for the glueball mass spectra come from computations on the lat- 
tice. Our backgrounds naturally give predictions for the = oo theory. Although 
there are large- A^^^ extrapolations (see for example |3^), there exist richer and more 
precise data for SU{3), especially for the excited glueball states. Therefore, we choose 
to fix our parameters in order to fit the available data for A^ = 3. We note that the 
error one makes for using A^ = 3 data instead of A^ = cxd is within 5 percent |37|. This 
is well within the error bars of the lattice computations for SU{3) (see pSi , ^). 



jPC 


Ref. I (m/y^) 


Ref. I (MeV) 


Ref. II (mro) 


Ref. II (MeV) 


Nc oo{m/^) 


0++ 


3.347(68) 


1475(30) (65) 


4.16(11)(4) 


1710(50)(80) 


3.37(15) 


0++* 


6.26(16) 


2755(70)(120) 


6.50(44) (7) 


2670(180)(130) 


6.43(50) 




7.65(23) 


3370(100) (150) 


NA 


NA 


NA 




9.06(49) 


3990(210)(180) 


NA 


NA 


NA 


2++ 


4.916(91) 


2150(30)(100) 


5.83(5)(6) 


2390(30)(120) 


4.93(30) 


2++* 


6.48(22) 


2880(100)(130) 


NA 


NA 


NA 


0-+ 


5.11(14) 


2250(60)(100) 


6.25(6)(6) 


2560(35)(120) 


NA 


0-+* 


7.66(35) 


3370(150)(150) 


NA 


NA 


NA 


R20 


1.46(5) 


1.46(5) 


1.40(5) 


1.40(5) 


1.46(11) 


Roo 


1.87(8) 


1.87(8) 


1.56(15) 


1.56(15) 


1.90(17) 


Rao 


1.52(8) 


1.52(8) 


1.50(5) 


1.50(5) 


NA 



Table 2: Available lattice data for the scalar and the tensor glueballs. Ref. I denotes |40| 
and Ref. II denotes 



I and [39|. The first error in the Ref.I and Ref. II correspond to the 
statistical error from the the continuum extrapolation. The second error in Ref.I is due to 
the uncertainty in the string tension \/a. (Note that this does not affect the mass ratios). 
The second error in the Ref. II is the estimated uncertainty from the anisotropy. In the 
last column we present the available large estimates according to [37|. The parenthesis 
in this column shows the total possible error following by the estimations in 



There exist a vast literature on the lattice computations for the glueball spectra. 

^30 listed the available data in 



We take as reference, the papers |38|, |39| and 



table ||. In that table Ref. I denotes and Ref. II denotes and |^ . Although 
we listed the lattice results also in the units of MeV, it is more convenient to use the 
units of ro (the "hadronic length scale") or ^/a (the confining string tension). In order 
to compare the data according to the two references, one should take oTq 1.36.^^ 



■^"^We thank H. B. Meyer, C. J Morningstar and M. Teper for pointing us to these references. 
■^^We thank H.B. Meyer for explaining this to us. 
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In order to avoid the error in the choice of the unit mass scale, we fit our param- 
eters by the mass ratios ratios, that we denote as: 



_ mo*++ _m2++ _mo-+ 

/too — ) ^20 — ; ^AO — • lO-J-J-J 

mo++ mo++ mo++ 



There is a shght mismatch for the values of these ratios in the refs. and [4C 



(see table I). Thus, in the next section, we shall present our results for fitting our 
parameters according to both of these references separately. 

Notice that we could have computed the string tension a by looking at the 
minimum value of the string frame scale factor, as explained in Section 3. To obtain 
any numerical information, however, would further require knowledge of the relation 
between the fundamental string tension and the AdS scale. The latter sets the 
overall mass unit. Since this relation is not fixed in our model it does not constitute 
an independent check. 

Finally we should mention that the experimental identification of glueballs in 
high energy experiments has a long and not very successful history. The main prob- 
lem is to find unambiguous criteria that would distinguish glueballs from others states 
(mesons, and hybrids) in the experimental data. Recent discussions on the status of 
the experimental glueballs search both for scalar and pseudoscalar ones can be found 
in references 

8.1.3 Fit for Reference I 
0++ and 2++ glueballs 



As we discussed above, the numerical integration of ( |6.5|) determines the spec- 
trum in terms of bo, up to a choice of scale. We showed that the mass ratios are 
independent of Aq and Ao. 

We fix Ao = 0.05, then vary bo to obtain the ratios i?oo = 1-87 and i?2o = 1-46 
(table I). We then fix the overall energy scale to set mo++ = 1475. As explained 
in Section 7, this is equivalent to fixing Aqcd, and completely determines the back- 
ground solution. We then compare our results with those in the third column of 
Table |. 

The value of 6o that fits Rqq = 1.87 is 6o = 4.2. Fixing this, we find R20 = 1.40. 
The masses for the lowest lying states are found to be: 

0++ mi, ma, ■ ■ ■ = 1475, 2753, 3561, 4253, 4860, 5416 ■ ■ ■ MeV. (8.12) 
2++ mi, ma, ■■■ = 2055, 2991, 3739, 4396, 5530, ■■■ Mel/. (8.13) 

Notice that the spectrum of excited states is in good agreement with the available 
data from Ref. I. We compare our results with the lattice data and the standard 
AdS/QCD predictions in Fig. ^. The glueball spectrum in the standard AdS/QCD 
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model is worked out in Appendix We note that the glueball spectra in the hard- 
wall model were first discussed in |^4[ and |^ with Dirichlet boundary conditions 
at the wall. 
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(a) 
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Figure 8: Comparison of giuebail spectra from our modei with bo = 4.2 (boxes), with the 
lattice QCD data from Ref. I (crosses) and the AdS/QCD computation (diamonds), for 
(a) 0^^ glueballs; (b) 2~^^ glueballs. The masses are in MeV, and the scale is normahzed 
to match the lowest 0^^ state from Ref. I. 
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Figure 9: The string frame scale factor in background I with bo = 4.2. The units on the 
horizontal axis are such that i = 1. 



String tension 



■^^There, we fixed tq by the meson data. If one leaves tq as a free parameter in the glueball 
sector, one can obtain better fits in the AdS/QCD set-up. For example, Q finds good fit with the 
Pomeron trajectory with Neumann boundary conditions. 
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/^From the first column of Table ^ we can estimate tlie fundamental string tension 
Tf in AdS units: 

Tff = a£2e-2^^(^*) = I!^±^e-'^^^''^ (8.14) 

The string frame scale factor is shown in Figure ^, and numerically we find that at 
the minimum e^"^"^^*^ ~ 2 x 10~^. This gives 

2^/ ^' = ^^6.24 ^ £^6.26 (8.15) 

The size of the UV geometry is several times the string length. This in particular 
shows that the dimensionless curvature invariant (in the Einstein frame) near the 
AdS^ boundary is 

i^R - -0.5 (8.16) 



Running coupling and QCD scale 

Using our choice of Aq and Aq we can compute the perturbative QCD scale as 
defined in eq. ( [7.10|) . We obtain 



Ap~ 290Me\/, (8.17) 

which is the correct order of magnitude, but not very close to the generally assumed 
value of around 200MeV^. This can be attributed to the fact that in deriving the 
latter value, the effect of five flavors of quarks is assumed to contribute to the running, 
whereas we are dealing with pure Yang-Mills. Moreover one should keep in mind that 
there is no unambiguous definition of A |^. It is more meaningful to compare with 
experiment the value of the strong coupling constant as at some energy scale. In 
order to do this, we must first identify the relation between our coupling constant A 
and the strong coupling constant a^. This is fixed once we set the parameter ^o, and 
can be obtained by comparing the one-loop beta-functions: 

r)\ 

-bo\^ + ..., (8.18) 



dhgE 



da. 00 9 ^11 / N 



dlogE 271 ' ' ""3 
^From the two expressions we find the relation: 
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where we have set A'^c = 3 in the last step^^. Using our numerical solution we can 
compute the value of as at some radius p; this can then be translated into an energy- 
scale through the relation E = EQe'^^'^\ in which £"0 is fixed by matching the lowest 
glueball mass (in our case Eq = 17630MeV^). For example, we find 



as{l.2GeV) = 0.34, 



which is within the error of the quoted experimental value |}47 
0.35 ±0.0134. 



(exp) 



(8.21) 
'1.2GeV) = 



+ glueballs: 

Having fixed Bq, we can now vary the parameter Ca to fit the lowest pseudoscalar 
glueball mass. First, we notice that for large values of Cq (greater than ~ 10), the 
spectrum depends very weekly on this parameter. This is shown in Figure p!0| 
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Figure 10: Lowest glueball mass in MeV as a function of c^. 

The ref. I value Rao = 1.52 is obtained for = 0.05. For this value the lowest 
pseudoscalar glueball masses are found to be: 



0" 



mi, ma, • • ■ = 2243, 3436, 4396, 4911, 5541, , • ■ ■ MeV. 



(8.22) 



Again, after fitting the lowest state, the first excited state is in good agreement 
with the lattice data of Ref. I. The comparison between our result and the data of 
ref. I is summarized in table ^ In figure |TT| we show the wave-function profiles of 
the lowest O"*"^, ^ and 2'^'^ states. 



■^^In Af — 4 SYM the identification is fixed by the D3 brane coupling to the dilaton, as = Qs ~ 

^■^The uncertainty on this value is not reported in p7| ; rather, it is an estimate obtained from 
the corresponding uncertainty in the data as{Mz) = 0.1202 ± 0.005 
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jPC 


Ref I (MeV) 


Our model (MeV) 


Mismatch 


iV, ^ oo IP^ 


Mismatch to — ^ oo 


0++ 


1475 C4%^ 


1475 





1475 





2++ 


2150 


2055 


4% 


2153 (10%) 


5% 


0-+ 


2250 (4%) 


2243 









0++* 


2755 (4%) 


2753 





2814 (12%) 


2% 


2++* 


2880 (5%) 


2991 


4% 






0-+* 


3370 (4%) 


3436 


2% 






Q++** 


3370 (4%) 


3561 


5% 






qH — h*** 


3990 (5%) 


4253 


6% 







Table 3: Comparison between the glueball spectra in Ref. I and in our model. The states 
we use as input in our fit are marked in boldface. The parenthesis in the lattice data 
indicate the percent accuracy. The data in |37| are given in terms of mass ratios; here they 
have been rescaled to match the lowest state mass. 



iA[r] 




r [mO] 



r[A] 40 



Figure 11: Normalized wave- function profiles for the ground states of the 0^^ (solid 
line),0 (dashed line), and 2^+ (dotted line) towers, as a function of the radial conformal 
coordinate. The vertical lines represent the position corresponding to E = mQ++ and 
E = Ap. 



8.1.4 Fit to Reference II 
0++ and 2++ glueballs 

We use the data from |^ for our fit, since that work includes values for some of the 
excited states masses. Although older, these data do not differ significantly from the 



more recent ones reported in where however no excited states masses are given. 
We fix bo to match i?oo = 1-54. The preferred value is now bo = 2.5. Then, we set 
the energy units so that mo++ = 1730. The lowest lying states have masses: 

0++ : mi, ma, ■ ■ ■ = 1730, 2697, 3321, 3853, 4319, 4747, 5139 , ■ ■ ■ MeV. (8.23) 



2++ : mi, ma, ■ ■ ■ = 2194, 2897, 3485, 3987, 4440, 4851, 5229, ■ ■ ■ MeV. (8.24) 
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Figure 12: Comparison of glueball spectra from our model with bo = 2.5 (boxes), with the 
lattice QCD data from Ref. II (crosses) and the AdS/QCD computation (diamonds), for 
(a) O"*"^ glueballs; (b) 2"*"^ glueballs. The masses are in MeV, and the scale is normalized 
to match the lowest 0"*"^ state from Ref. II. 



Running coupling and QCD scale 

Proceeding like in the previous section, we obtain: 

Ap = 356Me\/, as{l.2GeV) = 0.38 (8.25) 

These values are farther from the observational expectations compared to the result 
of the fit to Ref. 1. 



+ glueballs: 

We now vary the parameter Ca to fit the lowest pseudoscalar glueball mass. 





Ref 11 (MeV) 


Our model (MeV) 


Mismatch 


iV, ^ oo i7| 


Mismatch to Nc ^ oo 


0++ 


1710 (5%) 


1710 





1710 





2++ 


2390 (5%) 


2194 


8% 


2502 (10%) 


10% 


0-+ 


2560 (5%) 


2582 









0++* 


2670 (5%) 


2697 


1% 


3262 (12%) 


18% 


0-+* 


3640 (5%) 


3434 


6% 







Table 4: Comparison between the glueball spectra in Ref. II (taken from [^^) and in our 
model. The states we use as input in our fit are marked in boldface. The parenthesis in 
the lattice data indicate the percent accuracy. The data in [^] are given in terms of mass 
ratios; here they have been rescaled to match the lowest state mass. 

The ref. 11 value Rao = 1-50 is obtained for c = 0.023. For this value the first 
pseudoscalar glueball masses are found to be: 

0~+ mi, ms, ■ ■ ■ = 2597, 3434, 3965, 4407, 4851, 5246, • • • MeV. (8.26) 
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The comparison between our result and the data of ref. II is summarized in table 

i 



8.1.5 Dependence of the spectrum on the spectral parameter a 

Up to now we have set the spectral parameter a = 2, as it corresponds to linear 
confinement, oc n for large n. However, unlike in the case of mesons, there is 
no direct lattice or experimental evidence for such a behavior for the glueballs. In 
particular, the lattice simulations are only available up to n = 4 (for 0^"'^ only). 
Therefore, it is interesting to examine the dependence of the spectrum on a. We 
recall that the effective Schrodinger potential in the IR behaves as. 



V[r) 



r 



2{q-1) 



as r 



oo. 



(8.27) 



Hence, one expects that the mass spectrum will move upwards as one increases a. 
One also expects that the hard- wall approximation of AdS/QCD would correspond 
to a — i> oo. 
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Figure 13: The 0"'"'' spectra for varying values of a that are shown at the right end of 
the plot. The symbol * denotes the AdS/QCD result. 

We carried out the necessary numerical analysis for the 0^+ glueballs, for fixed 
values of Aq, &o and Aq and varying a. We fix bo = 4.2, as in the fit for Ref. I, 
so that the mass ratio of Rqo is 1.87 for a = 2. We normalize the spectra so that 
the lowest scalar glueball has the same mass for all a we consider. Our results are 

35 



depicted in fig.|13| where we also included the AdS/QCD result for comparison. 



^^To compare with AdS/QCD we fixed the value of Tq of [|j such that the first glueball lies at 
1475 MeV. 
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One indeed finds that as a increases the spectrum of our background approaches to 
that of standard AdS/QCD, and the agreement with Ref. I becomes worse for larger 
a. However, if we allow to change 60 we can fit the data equally well for a 7^ 2 but 
not too large, so there is no conclusive evidence that a = 2 is preferred. 

8.2 Background II: singularity at finite r 

In this section we compute the spectrum in a 5D background with different IR 
asymptotics, namely the one in which the IR singularity is at finite r. We assume a 
power-law IR singularity, 

y4(r) ~ (5 log(ro — r), r — Tq. (8.28) 

For the phase space variable, we take the same UV asymptotics whereas in the 

IR, according to eqs. ( |A.65|) , one has: 



X{X) = --^Q ^ . . . Q = ^v/rT5^ (8.29) 
As interpolating function we choose: 

To compute the spectrum we use the same procedure we employed in the previous 
example. We first integrate numerically the equations for the metric and dilaton, then 
we use a shooting method to find the mass eigenstates. We use 60 aiid 5 as fitting 
parameters. 

8.2.1 The glueball spectra in background II 

First, we obtain the spectrum for the same value of 60 that gives the best fit to the 
data from Ref. I, namely 60 = 4.2, and we vary the parameter 6?^ Varying 5 between 
5 = 1.01 and 5 = 10 we obtain the results in Table |^ a). To explore the dependence 
on 60 we fix 5 = 2 and vary 60 (see table |] b). For a wide range of 60 -Roo and i?2o 
are significantly smaller than the lattice values. 



0"*"^ and 2++ glueballs: Fit for Reference I 

To fit the data in Ref. I we use the following procedure: for different values of 5, 
we fix 60 to obtain the mass ratio i?oo = 1-87 as close as possible. Then we compare 
our finding for i?2o with the lattice value. Since the dependence on 60 for a given b 
does not follow a clear pattern, it is very hard to fit exactly any particular value of 
-Roo- It turns out that we were not able, with this ansatz for X(A), to obtain an i?oo 
larger than 1.65, for which i?20 = 1.3. 



36 



We always use b > \ because of the reasons discussed in Section 6.6.1 
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6 


Rqq 


i?20 




bo 


Rqq 


i?20 


1.01 


1.50 


1.20 




0.5 


1.47 


1.17 


1.05 


1.48 


1.19 




0.75 


1.42 


1.15 


1.1 


1.48 


1.19 




1 


1.39 


1.14 


1.5 


1.41 


1.16 




2 


1.38 


1.14 


2 


1.37 


1.13 




3 


1.37 


1.13 


3 


1.27 


1.09 




5 


1.37 


1.13 


4 


1.27 


1.08 




10 


1.37 


1.13 


5 


1.24 


1.07 




25 


1.40 


1.10 


7 


1.20 


1.05 




40 


1.41 


1.07 


10 


1.16 


1.04 




100 


1.47 


1.05 




(a) bo 


= 4.2 


(b) 


6 = 2 



= 1.87, R',o' = 1.46 



Table 5: Lowest glueball mass ratios for a) 6o = 4.2., /q = 0.05, for varying 6; h) 5 = 2., 
/q = 0.05, for varying bo 

0"*"+ and 2++ glueballs: Fit for Reference II 

Contrary to the case of Ref. I above, one can fit tlie value Roo = 1-56 in Ref. II 
(table H), by choosing bo = 0.96 and 6 = 1.01. However, we cannot find a set of 
parameters which also gives a good result for i?2o- For the aforementioned values of 
bo and 6, one obtains R20 = 1.25. 

8.3 Estimating the effect of the UV running 

In this subsection we investigate how the logarithmic running of the coupling in the 
UV affects the IR properties, such as the glueball mass spectrum. To address this 
issue, we compare the spectrum of background I (a = 2, 60 = 4.2) with another 
background obtained by keeping the same IR properties, but with a conformal fixed 
point in the UV. In the latter background, the geometry is asymptotically AdS^ up 
to power- law corrections, and the 't Hooft coupling fiows to a non-zero value A*, 
which can be chosen to be small. Such a geometry has the following asymptotics for 
the superpotential and /9-function in the UV (i.e. for A ~ A*): 

W,onf = Wo + Wi{X-X,f + ..., W^o = ^ (8.31) 
/?co„/(A) ~ -6oA.(A, - A) 60 > 0, A, < 1 (8.32) 
In the IR, we take the new background to have the same large A asymptotics as 



background I, as in (|8.2|) and (8^) with a = 2. Moreover, we fix the initial conditions 
and the parameter A* such that the strong-coupling scale of the two backgrounds are 
the same. As a definition of the strong coupling scale we take the slope of the scalar 
glueball mass spectrum: = A^n for large n. 



- 62 - 



As a simple example of an asymptotically conformal background with the desired 
IR properties we can take: 



e^(r) = -e (''Z^) , $(r) = $o + --^^ 1 + 3^ + - log . (8.33) 

One can easily check that the above solves Einstein's for an appropriate choice of 
superpotential, that is detailed in appendix 0. This is an example of an explicit 
"soft wall" background, in which both the metric and the dilaton are known exactly, 
and that can be derived as a consistent solution of Einstein's equations. 

We use the same shooting method as before to compute the mass eigenvalues. 
We can fix $o = $(0) and R in ( ^.33| ) to match the slope of the glueball masses 
found in the asymptotically free background. 

As an alternative background, we start with the exact superpotential: 

W^onf = W^o (l + \hl{\ - Kfyl^^ (9a + (26^ + 360 log [l + (A - X^)^""''' . 

(8.34) 

This amounts to a small modification of the superpotential ( ^.6|) , but it behaves 



asymptotically as (|8.31) in the UV 



The results are shown in figure |14 
8.4 Discussion 

Here we summarize the results of our numerical analysis. From the qualitative point 
of view, our general setup can reproduce the known features of the scalar and tensor 
glueball spectra. For example, as in the lattice studies, the O"*"*" states are lighter than 
the 2++ states, contrary to the AciS'/QCD models of P, §], in which the two towers 
are exactly degenerate. The pattern mi^^ < mi^"* seems to be a generic feature of the 
dual backgrounds in which the dilaton is taken to be non-trivial. We see numerically 
that this behavior is realized in all the backgrounds we considered, and it was also 



observed in [|T5]. Moreover, we always observe i?oo > -R20; which is common to all 
lattice results. 

The "linear" model with a = 2 seems to reproduce the pattern of excited spin-0 
glueballs found in the lattice study |^ which to our knowledge is the only work 



that computes the masses of such states. From the quantitative point of view, we 
can make the following comments. We remind the reader that our fits refer to mass 
ratios, as we can always choose arbitrarily the absolute energy scale. 

• For the infinite range background (background I) one can fit both sets of the 
available lattice data, Ref I. and Ref II, by fixing the parameter 60 • To check 
agreement with the lattice, one should look at the last column of Table 0, as 
our setup is supposed to describe 4D YM at large A^c- Notice that the large 



- 63 - 



40 - 



30 ^ 



20 - 



10 



■ ■ ★* 

■ ★ A ^ 



.It- 



■ ★ A ^ 
■ ★ A ^ 

■ ■*■ A 

■ ★ T ■A 

■ ^* A'^ 



10 



15 



20 



25 



30 



Figure 14: The comparison of the scalar glueball masses for the asymptotically free and 
the two conformal backgrounds: the stars correspond to the asymptotically free background 
(8.4) with 6o = 4.2 and Aq = 0.05; the squares correspond the results obtained in the 
background ( ^.33| ) with R = 11. 4i; the triangles denote the spectrum in the background 
given the superpotential (^.3l|) with bO = 4.2, Aq = 0.071 and A* = 0.01. These values are 
chosen so that the slopes coincide asymptotically for large n. 



Nc mass ratios i?oo and R20 are very close to the ones of Ref. I. Moreover, the 
uncertainties in Rqq and R20 for large Nc are larger than the ones reported for 
the glueball masses in both Refs. I and II. Our best fit for Ref. I is well within 
the large Nc error-bars. 

• The value of the spectral parameter a affects the results. We fix it to a = 2 in 
order to obtain a linear Regge trajectory. We note however that it is possible 
to fit the lattice data for a different set of values for bo and a ^ 2. In this case 
the large n asymptotics in the spectrum will not be linear. 

• As a general conclusion for the finite range background (background II), we 
can say that we could not find a range of parameters that yield good fits for 
both the scalar and tensor glueball masses. In particular, if one adjusts the 
parameters in order to fit the scalar ratio Rqq, then the tensor glueball masses 
turn out to be significantly lower than the lattice results, and outside the large 
Nc error bars. 

• We analyzed the dependence of the spectrum on the logarithmic running of 
the coupling in the UV, by comparing our results with a background where 
one has the same IR but a conformal fixed point in the UV. This background 
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has power law running for the couphng. One finds that for a fixed slope of 
the glueball spectrum, the overall scale of the masses do change. However it is 
possible to fit the lattice data by a choice of different parameters. Therefore, 
one can obtain in principle the same spectrum (at least for small n) in a theory 
where the UV is a conformal fixed point. 

• A final word on fitting the lattice data: our strategy is to fit i?oo by fixing the 
parameter 60 in our backgrounds and then obtain a prediction for the ratio 
i?2o- As we mentioned, this prediction falls into the error bars in the references 
I and II that account for the uncertainty in the large Nc limit (see table (^). 
Furthermore our predictions for the higher excited states also turn out within 
those error bars, if we assume the same large Nc uncertainty as for the lowest 
states^^ This is despite the fact that our method of fitting the data is somewhat 
crude. A better method would be to apply a global fit both for Rqq and i?2o- 
One expects from this method to produce better results for the higher excited 
states as well. 

Note added: while this paper was being completed we became aware of work 
along similar directions, which appeared in ||48[ 
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There is no large extrapolation available for the third and fourth excited 0++ states. 
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APPENDIX 



A. Characterization of confining backgrounds 

We consider the Einstein frame metric in the conformal coordinates, 

ds^ = e^^^") {dr^ + r]ijdx'dx^) , < r < tq, (A.l) 
where r = is the AdS boundary. The corresponding string frame metric is 

ds^ = e^^^^'-) {dr^ + riijdx'dx') , Asir) = A{r) + '^^r). (A.2) 

Given the behavior of the scale factor close to the singularity, the asymptotic behavior 
of the dilaton is uniquely fixed by the first of eqs. 



^\r) = -^(^A{r)-A\r)). (A.3) 

Knowledge of A{r) and $(r) uniquely determines the asymptotics of the phase space 
variable X, therefore those of the /3-function from eqs. (P^TBf). Asymptotics of the su- 



perpotential W can be determined from eqs. ( p. 18 ), or from the second eq. in (|2.12| ). 



X, (3 and W can then be expressed as functions of $ by inverting asymptotically the 
relation between $ and r^®. 

Therefore, we can parametrize different backgrounds by the asymptotics of the 
scale factor alone, since this completely determines the asymptotics of all other quan- 
tities. The singularity can be at a finite or an infinite value in the conformal coor- 
dinate. We discuss these two cases separately. For all cases analyzed below, we give 
the IR asymptotics of the following quantities, found by the following equations: 

• Einstein frame scale factor A(r), 

• Dilaton and 't Hooft couphng $ = log A. 

• String frame scale factor : 

As = A+'^^ (A.4) 

• Einstein frame and string frame curvatures^^ : 

R ~ e-2^i^ Rs ~ e-2^M| (A.5) 



■^^This can be done in backgrounds where the NEC is satisfied, see Section 2. 

■^^In the Einstein frame there are two independent curvature invariants, {dr^)^ and the Ricci 
scalar. They both behave asymptotically as e~'^^A^, and will be denoted collectively by R. The 
same holds for the string frame. 
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Phase space variable and /3-function, 



X(<|.) = ^, /?(A) = 3AX(A) (A.6) 



Superpotential 



W{X) ~ exp 



(A.7) 



Dilaton potential (in Einstein and string frame actions): 

vm = -I (^) ' + vsm = e-^*/v($) (a.s) 

Metric and dilaton asymptotics in the domain-wall coordinate u: 

u = I dre^ir) (A.9) 



A.l Unbounded conformal coordinate 

If the space-time extends over an infinite range of the r coordinate, the Einstein frame 
scale factor e'^(r) necessarily vanishes as r — ^ oo, as a consequence of eq. ( |2.24| ). 



Therefore, A{r) — oo as r — oo. We analyze two possible types of behavior for 
A{r), logarithmic and power law (the latter was also discussed in Section]^). In both 
cases the singularity is at a finite value uq in "domain wall" coordinates. 

A. 1.1 Logarithmic divergence 

Consider backgrounds such that, for large r: 

A(r) ~ -7logr + . . . 7 > 1. (A.IO) 

The constraint 7 > 1 comes from the Null Energy Condition discussed in Section 2. 
7 = 1 corresponds to AdS asymptotics in the IR, which does not lead to confine- 
ment. For any 7 > 1, there is no confinement either, as we show below. We have, 
asymptotically: 

(A.ll) 



From ( |A.3|) we obtain: 



Eq. ( |A.12| ) integrates to: 



$^ ~ (A.12) 



2 

$ ~ - A/72 - 7 log r. (A. 13) 
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iFvom eq. (|A.2|) , the string frame scale factor behaves asjTiiptotically as: 

Asir) ~ -7 (l - v^T^T^) logr, (A.14) 

Since the overall coefficient is negative (7 > 1), As{r) —00 as r —>■ 00. Therefore 
the string tension vanishes and there is no area law in this case. These are the 
asymptotics of the relevant quantities: 

as r — > 00 : 



A ~ -7 logr, 7>1; Q = ^^l_l<^ (A.15) 

color confinement: NO (A. 16) 

3 9 

^ -^y^^ - ^\ogr = -'jQlogr, (A.17) 

Asr^—f(l-lQ]logr, (A. 18) 



2 

~ r2(T-i) ^ cx), (A. 19) 

R, ^ ,2(.-i)-3.g f 1 < 7 < i(l + v^) (A.20) 

\oo 7>|(1 + V5) ^ 

X{\)r^-^^, Wi\)r^\Q (A.21) 

V ~ A^^, Vs = X'W ~ A^^-t (A.22) 

ur-^uo-O (-1^) (A.23) 

A(u) 2_log(Mo-M). (A.24) 

7 - 1 

A. 1.2 Power-law divergence 

Next we consider the following large r behavior: 

A(r) ~ -Cr" + ..., C>0,a>0, (A.25) 

where the precise nature of the subleading terms is immaterial. This case was dis- 
cussed in Section (^. It leads to confinement if and only if a > 1. We have: 

A ~ -Car"-\ A ~ ~Ca{a - l)r"-2 (A.26) 

Notice that A/A ~ r^^, therefore eq. ( |A.3| ) is solved, asymptotically, by: 

<^ = -lA+hog\A\ + % + o(-] (A.27) 
2 4 \r J 

where we have kept the first subleading term, which is universal and independent of 
the subleading terms in (|A.25|) . The string frame metric, from eq. ( |A.2D , is: 

As ~ ^^y^ logr + ho + (A.28) 
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Notice that the leading terms cancel. ( |A.2S ). Therefore: 

{-oo, < « < 1 
const, a = l (A.29) 
+00, a > 1 

and we have confinement if and only if a > 1. In the borderline case a = 1, As 
asymptotes to a finite constant as r 00. The string frame metric is asymptotically 
Minkowski, and the dilaton is linear in r, up to subleading corrections. The string 
of minimal world-sheet area stretches all the way to r = 00, but the confining string 
tension is nevertheless finite. 

We list below various relevant quantities: 

as r — >• 00 : 

A ~ -Cr", a > 0, C > 0; P = < 1 (A.30) 

a 

color confinement: if a > 1 (A.31) 

$~^Cr" + ^(a-l)logr, (A.32) 

^5~^^^logr, (A.33) 
R ^ e^^'-V^^"-^) ^ 00, (A.34) 

(A.35) 

X(A) ~ (1 + ^1^) ' ^(^) ~ (logA)-^Ai (A.36) 
~ (logA)^A3, Vs = X'W ^ {\ogXf (A.37) 

The domain wall coordinate u terminates at a finite value uq, as the integral in eq. 
( |A.9|) converges as r — > 00. The metric and dilaton in this frame are, close to the 
singularity: 

u Mo, log(Mo -u) Cr" (A.38) 

A{u) ~ log(Mo -u) + P log [- log(Mo -«)] + ..., (A.39) 

3 3 
~ -- log(uo -u)--P log [- log(Mo - u)] (A.40) 

A. 2 Finite range of the conformal coordinate 

Now suppose that the singularity is at a finite value of the conformal coordinate, 
r = tq. By monotonicity of A{r), the scale factor at the singularity either vanishes, 
or stays finite. 

A. 2.1 Finite A{ro) 

If v4(ro) is finite, the singularity must be caused by non-analyticity in A. The dilaton 
may stay finite at vq, or asymptote to +00 (we are assuming strong coupling in the 
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IR, so we exclude the case $(ro) = —oo). In any case, the string frame scale factor, 
A + 2$/3, is either finite at ro or asymptotes to +00, therefore it must have a 
minimum for some r* in the range (0, Tq]. The value at the minimum must be finite 
(otherwise there would be a singularity at < tq) , leading to a confining string with 
non-zero tension. 



According to the identification (|2.15|) , the fact that the Einstein frame scale 



factor is nowhere vanishing means that the dual 4D theory is defined only above a 
certain energy Emin ~ e"^""". We will discard this case for a different reason: there 
is no screening of the magnetic color charge. 



A. 2. 2 Power-law divergence 

Next, we consider the case when the Einstein metric scale factor vanishes at some 
r = ro as a power-law: 



C 

A(r) a > 0, C > 0. (A.41) 

Below we show that the string has a finite tension for all acceptable values of a and 
C. The argument we present holds for any generic subleading behavior. One can 
easily check that the solution of (|A.3|) close to ro is given by 



$(r) ~ -^A(r) + ^ log |i(r)| + $0. (A.42) 



This ansatz solves (|A.3|) up to a term proportional to [A/Ay ~ (ro — r) ^, which 



for a > is subleading w.r.t the term A^ ~ (ro — r)^""*"^ in eq. ( |A.3| ) . The string 
frame metric asymptotes as: 



/Ig ~ ^ log i ^ - ^ log(ro - r). (A.43) 



The leading terms cancel, and the first subleading term is universal. Eq. ( |A.43| ) 
shows that Ag +00 as r ^ ro for any positive a, and we always obtain a confining 
string. 
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We list below various relevant quantities: 



as r ro : 

, a>0,C>0; P = ^^^>1 (A.44) 

color confinement = YES (A. 45) 

$ ~ ^^^r^ - + 1) log(ro - r), (A.46) 

^5~-^^^log(ro-r), (A.47) 
1 



R ~ 7 T^T^eC-o-)" ^ oo, (A.48) 

(ro - r)2("+i) ^ ^ 

(cx) < a < 1 

const a = l (A.49) 
a > 1 

X(A) ~ (l + ^i^) ' ^(^) ~ aogA)^ At (A.50) 

\/~(logA)^Ai l/s = A-ty ~ (logA)^, (A.51) 

~ log(uo -u) + P log [- log(Mo - m)] + . . . , (A.52) 

3 3 

~ -- log(Mo -u)--P log [- log(Mo - u)] . (A.53) 

A. 2.3 Logarithmic divergence 

In this case we have, asymptotically: 

A ~ 51og(ro -r), 5>0, (A.54) 

and 

5 S 

-, Ar^— -. (A.55) 

(^0 - r) (ro - ry 

From ( |A.3| ) we obtain: 

* - 4 (^^' (^-^^^ 

Eq. ( |A.56| ) integrates to: 



$ V52 + 51og(ro - r). (A.57) 

where we chose the branch ($ > 0). The string frame scale factor behaves asymp- 
totically as: 

As{r) ~ 5 (^1 - + log(ro - r). (A.58) 
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For large r it asymptotes to +00, as the overall coefficient is negative for positive 6. 
Thus, the fundamental string confines. 
In this case we have: 



as r — i> To : 



2 / 12 

A~51og(ro-r), 6>0; Q = ^y^ + -^>^ 

color confinement = YES 

9 



$ ~ -VW+6\ogr = -6Q\ogr, 

As ^ —S ( 1 — -Q] logr = — — logfrn — r) 

* V ^ 1 + 3Q/2 " ' 

^ ^ _ ^)2(<5+l) ^ ^' 

3Q 

-Rs ~ ('^o - r) 3Q/2+1 ^ 00 



M~Mo-0((ro-r)^+i) 

+ 1 



fA 


59) 


(A 


60) 

\J\J J 


(A 


01) 


(A 


62) 


(A 


63) 


(A 


64) 


(A 


65) 


(A 


66) 


(A 


67) 


(A 


68) 



B. Magnetic charge screening: the finite range 

Here we want to determine the potential between two magnetic charges at large 
separation, for the type of backgrounds with ro < +00. The case ro = +00 was 



treated in section (|3.4j) . 
B.l A{ro) finite 

If the Einstein frame scale factor does not vanish at the IR singularity, the D-string 
frame scale factor cannot vanish either, and there is no difference between the cal- 
culation of the 't Hooft loop and that of the Wilson loop on the same background. 
As explained in Appendix (^), Section A. 2.1, in this case the electric string confines. 
Therefore the magnetic string confines too. These kinds of background fail to satisfy 
an important test for a candidate holographic dual of QCD. The same consideration 
applies to all theories where the 5th dimension terminates at a regular IR boundary. 

B.2 A{ro) -00 

We treat the case of power-law decay of the scale factor e^. The exponential case 
can be discussed along the same lines. We take 

A ~ 51og(ro -r), 5 > 0. (B.l) 
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From eqs. (|A.54 ) and (|A.57 ) we see that the D-string scale factor is asymptotically 
(as r ~ ro) 

An = A+^r^5(^l- + log(ro - r). (B.2) 

For 6 < 1/15, — » +00, the scale factor diverges at the singularity and the 
magnetic string confines. For S > 1/15 the scale factor vanishes as a power-law: 



e'^^« ~ (ro - r)^, j = 25 (^1 - ^VTTS^^ > . (B.3) 

In this case the magnetic string tension vanishes. To investigate the potential between 
two monopoles at large L, it is sufficient to translate our setup into the notation of 
[0 and use their results: defining s = ro — r, we are in the situation described in 

0, with /(s) = g{s) ~ s"^ as s — 0. In their notation, this is the case /(O) = and 
k =< j + 1 (since k = j = 'y). From their general analysis it follows that, for small 
sp (the turning point of the world-sheet), 

L{sf) ~ 4 (B.4) 

1. e. L{sf) vanishes as sp approaches the singularity. The same is true in the UV: 
L{sp) always vanishes close to an asymptotic AdS region. Therefore, it must be that 
L has a maximum value Lmax for some Vmax < ^o, and there is no smooth solution 
of the geodesic equation for L > Lmax- As we argued earlier in the case of infinite 
ro, the magnetic charges are free for L > L^ax- 

The behavior of L{rp) in the case with exponential fall-off close to ro cannot be 
deduced directly from the results of IT^, but it can be addressed by adapting the 



discussion in Section 3.4, and the result is the same, i.e. L{rp) cannot diverge. 



C. Fundamental string world-sheet embeddings in the pres- 
ence of a non-trivial dilaton 

The relevant world-sheet action is 

^ = ^ 1 rf'e ^g''^G,,{X)d.X>^dpX'' + d'^^R^'^^X) , (C.l) 

Instead of solving the equations we will do a simpler test. We will show that the 
contribution of the dilaton coupling to the full energy of the string is negligible in 
the limit where the distance between the endpoints of the string becomes large. 

We evaluate the action in the vicinity of the point r = r^, at which the scale 
factor of the target space metric has a minimum. We use the conformal coordinate 
system: 

cis2 = e2^=M(cix^ + rfr2), (C.2) 
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where As is the string frame scale factor, As = A + |$. We assume that the 
contribution of the second term in (|C.1|) is small with respect to the first term and 
we confirm our assumption, a posteriori. Then the leading term in the solution to 
the equation of motion that follows from ( C.l ) is. 



gab = gab = G^.daX^d.x". (C.3) 

We fix the diffeomorphism invariance on the world-sheet by choosing r = X", a = Y . 
Here, Y is the direction in the Minkowski space on which the quark pair lies. Using 
( |C.2| ) and ( |C.3| ), it is straightforward to compute the Ricci scalar on the world-sheet. 
One finds, 

^gR = ((1 + B')B'A': + BB'A'^) , (C.4) 

(1 + 5(r)2)2 

where we defined B{r) = {dy{r)/dr)~^ and the primes denote derivatives w.r.t. r. 
Notice that i? = at the world-sheet turning point. The second term in ( |C.1| ) 
becomes, 

^(2) =T [ \ ((1 + B')B'A', + BB'A:) $(r). (C.5) 

We assume that the scale factor Ag has a minimum at a point r,,. When the world- 
sheet turning point reaches tq. A' and B in (|C.5|) both vanish and A" and $ are some 
positive constants and the quark pair distance L = J dy diverges. Then it is clear 
from above that, 

S{2) const, 



therefore it is bounded in L, whereas the Polyakov term in (|C1D diverges linearly in 



L (under the aforementioned assumptions). Hence we can ignore the dilaton coupling 
in ( |(J.1| ) consistently. However, one has to be careful about the situations in which 
the integrand in ( |(J.5| ) asymptotes to a constant. In these cases, S(2) oc L and one 
cannot ignore the dilaton corrections to the induced metric. 

The picture we assume is as follows: the string world-sheet is smooth with a 
single turning point at rt. The geometry of the string is determined by a single 
boundary condition that we can take as the length between the end-points of the 
string on the boundary, L. As L is made larger the turning point approaches 
the minimum of Ag that we call r^. In particular we are assuming that there is a 
single minimum for Ag. As L approaches infinity, the greater part of the world-sheet 
falls into the minimum r^,. This picture is valid for all of the backgrounds that we 
analyzed in this paper. Thus, indeed the only term that has a potential divergence 
is the first term in ( |(J.1|) . 
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D. Singularities of the tachyon 



Here we analyse the properties of eq. ( |6.63| ) in the case with the following asymptotics 

A,(r)~^ log r/i?, «>1 (D.l) 

First, assume r(r) is nonsingular for any finite r. We want to analyse the be- 
havior near r = oo. Asymptotically, ( |6.63| ) becomes: 

We are interested in solutions that diverge as r — cx). First, suppose that r — » 00, 
but f stays finite. In this case, the third term in eq ( p.2| ) would be much larger than 
all others, and the equation would not be solved asymptotically. Then we conclude 
that as r ^ 00, f ^ 00 as well. Then the last two terms dominate eq. ( p.2| ), and 
the solution behaves as: 



r(r) ~ To exp 



2R 

a&^ 



00. (D.3) 



where tq is an integration constant. 

Now we want to check if it is possible for the tachyon to diverge at some finite 
value r=K, where the metric and the dilaton are non-singular. Then, close to r* : 

As = Ao + {r- r,)A, + ^(r - r,)^ + . . . (D.4) 

$ = $o + $i(r-n) + ^(r-r*)2 + ... (D.5) 

and we can approximate eq. ( |ti.63D by: 

f + (3Ai - $1) f + e^^V'^ + e-'^" [4Ai - $1] [if + /iV(f)2 = 0. (D.6) 

If r — >■ CX3 at r*, then the ratios f/r, f/f and f/f all diverge at r*, implying that 
the terms in eq. ( p.6| ) proportional to f and (f)^ diverge faster than all other terms. 
Therefore close to we can further approximate eq. ( p. 61 ) by: 

r + e-2^« [4Ai - $1] {jf = 0. (D.7) 

This equation is solved by: 



1 e 



2Ao 



which is not consistent with the assumption that r diverges at r^,. Notice however 
that the approximation we made in writing eq. ( p.7| ) still holds if we make the 
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weaker assumption that f, and not necessarily r, diverges at r*. Then, eq. (p.^ ) 
correctly describes the asymptotics near r*. This is, in fact, the generic behavior for 
arbitrary boundary conditions, the point r^, being determined by initial conditions. 

There is one situation when the above argument breaks down, i.e. when there 
exists a point r^,^, at which AAi — $i = 0. In this case the term in ( |D.6D proportional 
to {fY acquires an extra (r — r**) factor, and it is possible to solve the equation 
asymptotically with the last two terms: 

(r-r„)e-2^« [4^2 - $2] iff + fiMif ~ ^ r ~ (r-r„)l+^ h = 

(D.9) 

If 1 + /i < 0, this is consistent with r(r) diverging at r^,*. 



E. The superpotential versus the potential 

In this appendix we analyze the (important) details of passing from the potential to 
the superpotential. The first order non-linear differential equation that relates the 
two is 



For a given choice of V^(A) there is in principle a one-parameter family of solutions 
to (|E.1|) . This extra parameter that corresponds to the initial condition of the first 
order differential equation is compensating from the fact that the equations of motion, 
written in terms of the superpotential are first order and therefore have one less initial 
condition. 

According to our previous discussion, the different solutions to ( [E.l| ) correspond 
to different /5-functions, and hence to different 4D gauge theories. However, here 
we show that the requirement of color confinement in the IR singles out a unique 
superpotential among these. 

As described in detail in [0] asymptotic freedom implies that for small A: 

00 

V = Y,yny' = vq + Vi\ + o{\^). (e.2) 

n=0 

A straightforward calculation shows that the one-parameter family of solutions to 
the non-linear equation ([E.l|) has the structure 

^ n=0 ^ n=0 ^ n=0 

where C 7^ is the one free parameter, while all other coefficients are uniquely fixed 
by the expansion of the potential: 

27 27 27 27 

Co = , Ci = Vi , C2 = V2 , C3 = 1/3 (E.4) 

° 256 ' 352 ' ^ 448 ' ^ 544 ^ ' 
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27 27 27 27 27 

^^= 640^^ ' ^^= 736^^ ' ^^= 832^^ ' ^^= 928^^ ' ^^= 1024^^ 

" 19 ' 176 ' 50 ^ ' 

^ 42CI + 75C1C3 + 48C0C4 ^ 57C2C3 + 48C1C4 + 3OC0C5 

= 112 ' = 62 ^^-^^ 

E,^^C^,C. , ^^^ 3421CoC,- + 2128CgC. 
" 171 ° ' 16720 ^ ' 

There exists however a second branch of solutions again parametrized by a single 
integration constant C, disconnected from the family Wc{)\)'- 



WciX) = Wo + W^X + OiX^) + CXfe-^- + ■ ■ ■ , Wo = sJ^Vo, = y , 

(E.9) 

where the dots denote terms that vanish faster than exp(— 4/6oA). ( [E.3| ) together 
with ( [E.9|) provide the general solution of ( [E.l|) . All the solutions of the continuous 
family are singular for A ^ 0, and the corresponding metric is not asymptotically 
AdS^. AdS asymptotics require that W const in the UV and therefore single out 
the branch. 

It is the one that leads to the correct /3-function, at least perturbatively. However 
one wonders how the value of C is fixed among the solutions of this branch. For this 
purpose we turn to the analysis of the IR asymptotics below. 

We investigate the same problem, namely the solution of equation ( [h].l[ ) but in 
the IR region. We have argued that confinement requires that the IR asymptotics of 
the superpotential must be of the form, 

V ^ A^^e^Q*, $ ^ 00. (E.IO) 

with 1 < Q < 4/3. Solving the equation for the superpotential we find the one- 
parameter family of solutions 



Wc = (i) ' {Ce^'" + ^(4 - 3g)-^e(2«-t)*<|.^| + ■ ■ ■ , (E.ll) 

that is parametrized by one integration constant C. In addition there is the "singu- 
lar" solution. 

Consider now the IR asymptotics of the potential that we expect to give confinement 
in the case where the range of the radial coordinate is infinite. 

V = Ael'^ (E.13) 
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From ( [E.ll ) and ( Pj^.l2|) it is the singular solution ~ es* ^-^/^ that we need, rather 



than the continuous family. 

Now, one can show by scanning the space of solutions with all possible asymptotic 
behaviors in the IR, that the above leads in the UV to W^^ with C = 0. This is 
done by solving ( [li).l| ) numerically. We observe that the requirement of confinement 
singles-out the unique singular solution out of an infinite set. 

In view of this phenomenon, the practical strategy for choosing the correct po- 
tential is to first choose the appropriate superpotential (or the function X) and then 
use the differential equation (ETI) to compute the relevant potential. 



F. Standard AdS/QCD Glueball spectrum 

In this appendix we consider the standard AdS/QCD model where the background 
geometry is AdS^ with an IR cut-off at r = tq. The dilaton is constant. The glueball 
spectra in this model were first computed by employing Dirichlet boundary conditions 



in the papers and |^5|. However, here we shall allow for more general boundary 
conditions. 

In this geometry, both the scalar and spin-two glueballs spectra are determined 
by the following equation: 

^--i + m'^ = 0, (F.l) 

r 

The corresponding effective Schrodinger potential is. 



Vs = r < ro (F.2) 



and there is an infinite wall at r = tq. 

The solution to ( [F.l| ) that is normalizable in the UV is, 

^ = r^J2{kr) ^r^ as r^O. (F.3) 

The important difference between our backgrounds and AdS/QCD is that in 
AdS/QCD the normalizability condition in the IR does not restrict the spectrum. 
Indeed all the solutions of (|F?1|) with the UV asymptotics ( [I''.3[ ), are normalizable 



in the IR. What discretizes the spectrum is the boundary condition at r = tq. In 
general this can be a mixed boundary condition that may be written as, 

i{ro) - Q e(ro) = (F.4) 

Here Ct are real numbers and one can have different Ci for different particle species 
i. Therefore, the free parameters to fit the data are ro, Co++ and C2++. In the 
standard AdS/QCD model, the value of ro is determined by fitting the pion mass 
which yields ro = 1/322 MeV'K 
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We want to determine Co++ and C2++ to obtain a best fit to the lattice data. 
To obtain a best fit to the first O"*"^ glueball (1730 MeV), one has to avoid the first 
solution that is shown in fig. We note that linim^o ^/^("^^o) = 2. Hence one 
needs Co++ > 2. A quick glance at the fig.^ shows that the best fit (highest possible 
mass) for the first 0+"*" mass is obtained by setting Co++ = 2 + e in the limit e — 0+. 
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Figure 15: Plot of ^/^ as a function of mrQ. The spectrum is determined by the points 
that correspond to the intersection of this plot and the horizontal line ^/^ = const. 

Then one determines the O"*""*" masses as, 

mi, m2, ■ ■ ■ = 1651, 2710, 3734, 4764 5778, 6792 ■ ■ ■ MeV. (F.5) 

However now the best fit for the 2+"*" masses is given by the same IR boundary 
condition in the IR, i.e. C2++ = Co++ This gives the same mass spectrum for the 
spin-2 gluebaUs as in ( [b'.SI ). 

G. A simple analytic solution with AdS and confining asymp- 
totics 

We present here a simple analytic solution for the metric and the dilaton that inter- 
polates between AdS^ and a constant dilaton near the UV boundary and an Gaussian 
confining behavior and strong coupling in the IR. It is given by 

e^(r) = -e-(^/«) , $(r) = $0 + - — ^ 1 + 3— + - log j= . (G.l) 

''^^There is of course the possibility of choosing Co++ bigger than C2++ > 2. Ho'wever then the 
first scalar glueball masses becomes smaller than 1651 Mev and the discrepancy with the lattice 
data increases. 
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There are two dimensionless parameters that characterize this solution. The 
first is the asymptotic value of dilaton $o in the UV parametrizing the UV gauge 
coupling. To mimic YM, it should be taken large and negative. The other is the 
ratio -R/£, that characterizes the confinement scale R. 

The associated superpotential can be given in implicit form (as a function of r): 

W{r) = ^e^''/""' (1 - 2rVi?2) (G.2) 

Derivatives with respect to $ can be computed using the chain rule, W = W/^ etc. 
Note that this solution may look similar but is very different from the so called soft- 
wall AdS/QCD model. Although the dilaton has the same behavior, the metric here 
in the IR is very different from AdS. On the other hand the metric in the soft- wall 
AdS/QCD model is AdS everywhere. 
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